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Abstract 



Central to the development of any new theory is the investigation of the observ- 
able consequences of the theory. In the search for quantum gravity research in phe- 
nomenology has been dominated by models violating Lorentz invariance - despite 
there being, at present, no evidence that Lorentz invariance is violated. Causal set 
theory is a Lorentz invariant candidate theory of quantum gravity that seeks not to 
quantise gravity as such, but rather to develop a new understanding of the universe 
from which both general relativity and quantum mechanics could arise separately. 
The key hypothesis is that spacetime is a discrete partial order: a set of spacetime 
events where the partial ordering is the physical causal ordering between the events. 
This thesis investigates Lorentz invariant quantum gravity phenomenology motivated 
by the causal set approach. 

Massive particles propagating in a discrete spacetime will experience diffusion in 
both position and momentum in proper time. This thesis considers this idea in more 
depth, providing a rigorous derivation of the diffusion equation in terms of observable 
cosmic time. The diffusion behaviour does not depend on any particular underlying 
particle model. Simulations of three different models are conducted, revealing be- 
haviour that matches the diffusion equation despite limitations on the size of causal 
set simulated. 

The effect of spacetime discreteness on the behaviour of massless particles is also 
investigated. Diffusion equations in both affine time and cosmic time are derived, and 
it is found that massless particles undergo diffusion and drift in energy. Constraints 
are placed on the magnitudes of the drift and diffusion parameters by considering 
the blackbody nature of the cosmic microwave background. Spacetime discreteness 
also has a potentially observable effect on photon polarisation. For linearly polarised 
photons, underlying discreteness is found to cause a rotation in polarisation angle and 
a suppression in overall polarisation. 
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Notation 



-< The 'precedes' relation on a causal set. 

-<* The link' relation on a causal set. 
d(x, y) The length of the longest chain between elements x and y on 
a causal set. 

x The derivative of x with respect to time. 

A*, A\ *A The complex conjugate, adjoint, and Hodge dual of A, 
respectively. 

M 4 4-dimensional Minkowski spacetime with metric rj^ u . 
H 3 The mass shell, momentum state space for massive particles. 

Hq The massless particle momentum state space. 

B The Bloch sphere polarisation state space. 

A metric signature ( — h ++) is used throughout. 

Planck units with c = h = G = 1 and Boltzmann constant ks = 1 are used. 
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Chapter 1 



Introduction 



... I shall admit a system as empirical or scientific only if it is capable of 
being tested by experience. 

Karl Popper 

Is there really a singularity inside a black hole? Did the universe have a beginning? 
These are just two of the many questions that can't be answered until a theory of 
quantum gravity is developed. Quantum gravity isn't demanded by any unexplained 
experimental results. The predictions of general relativity and quantum theory agree 
with experiments wherever they have been tested. The main motivation behind the 
long standing search for quantum gravity is a philosophical desire for a single unified 
approach to physics. General relativity and quantum theory provide entirely different 
world views: is the universe a dynamic four dimensional Lorentzian manifold, or 
is there a fixed background spacetime on which fields exist? Or is reality stranger 
than either of these options? Not only are our main theories of physics seemingly 
incompatible, they contain within themselves the signs of their own breakdown. The 
singularity inside a black hole is not so much a prediction of general relativity as an 
indication that we are attempting to use the theory in a realm where it doesn't apply. 
Renormalisation and the associated problems in quantum theory suggest to me, at 
least, that a key part of our understanding is lacking. 

The search for quantum gravity has become not merely a search for a mathemat- 
ical method of tying the theories together, but a search for a completely new under- 
standing of the fundamental nature of spacetime. A number of approaches to quan- 
tum gravity exist. Most, including the prominent approaches of string theory and 
loop quantum gravity, give some preference to quantum theory, attempting to liter- 
ally 'quantise gravity' to fit it into a unified framework. Others, such as the causal 
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set approach that will be discussed in this thesis, begin with a completely new fun- 
damental structure and hope that general relativity and quantum theory will arise in 
the appropriate limits of the theory. No approach to quantum gravity as yet, can be 
considered complete. 

This thesis will focus on the possibility that spacetime may be fundamentally dis- 
crete and, more importantly, whether such a hypothesis is testable. If the explanation 
of experimental results has not yet required any theory of quantum gravity how are 
we to test the many theories that are being developed? If no effort is made to test 
developing theories they risk becoming only interesting mathematical constructs and 
losing any connection with a description of reality. As research areas expand and it 
becomes more and more difficult for an individual to comprehend an entire subject, it 
is crucial to remind ourselves that if we claim to be physicists we should be attempting 
to check our theories against the real world at every step. 

Since theories of quantum gravity must reproduce the results of general relativity 
and quantum theory in all areas where they have been shown to hold, and introduce 
new phenomena only on very small spacetime scales (or large energy scales), directly 
testing them is, at least in the near future, impossible. To test quantum gravity we 
need to seek ways in which such small effects could become amplified and result in 
deviations from the standard predictions of general relativity and quantum theory. 

It may seem strange to seek a signature for fundamental spacetime discreteness in 
large scale phenomena. Consider, however, the discreteness of matter. When Einstein 
provided an explanation for Brownian motion, he provided the last piece of evidence 
necessary to convince any doubters of the atomic nature of matter. That matter was 
made of atoms and molecules was, of course, not a new idea at that point. Einstein did 
not resolve the question of the discreteness of matter by directly observing atoms or 
molecules, but by recognising that the already observable phenomenon of Brownian 
motion demonstrated their existence. Likewise, to determine if spacetime is discrete 
we should not look to magical future technology that may let us see spacetime 'atoms' 
directly, but rather seek some currently observable phenomena that reveals the an- 
swer. 

Investigations into quantum gravity phenomenology have thus far focused over- 
whelmingly on producing and observing small violations of Lorentz invariance, often 
by introducing modified dispersion relations. As yet there is no evidence that Lorentz 
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invariance is violated at all. This is not to say that Lorentz violating quantum grav- 
ity phenomenology should not be explored: Lorentz violations can be investigated 
by current experiments and thus offer a very useful way of testing quantum gravity 
theories. It is often unclear, however, whether Lorentz violations are a necessary con- 
sequence of any particular quantum gravity theory, or simply a possible outcome that 
happens to be the only one that can currently be tested. Observations have forced very 
tight constraints on Lorentz violating models, for a recent review see 0. Violations 
of Lorentz invariance can never be ruled out, as experiments can only constrain the 
effects to be smaller. While it is important to explore these constraints to experimen- 
tal limits, it is clear that Lorentz invariant quantum gravity phenomenology deserves 
more attention than it currently receives. This thesis attempts to begin to remedy this 
problem. 

Causal set theory, a Lorentz invariant, discrete theory of quantum gravity, provides 
the primary motivation for the phenomenology discussed in this thesis. Although the 
work in this thesis is discussed within the framework of causal set theory, it should be 
noted that many of the conclusions are more generally applicable. 

To investigate Lorentz invariant phenomenology, I focus on the effect spacetime 
discreteness would have on the propagation of particles. Particles travelling through 
a spacetime with an underlying discreteness are expected to deviate from the contin- 
uum geodesies predicted by general relativity. For massive particles this was first con- 
sidered by Dowker et al. |3]], who found the propagation of massive particles could 
be described by a diffusion equation in the proper time of the particle. This thesis 
extends the work of Dowker et al. to derive a diffusion equation in observable lab- 
oratory time ('cosmic time') and also numerically investigates causal set models for 
massive particle propagation. A phenomenological model for massless particles trav- 
elling in discrete spacetime is also developed and astrophysical and cosmological data 
allow the free parameters in the model to be constrained. The models discussed here 
do not, unfortunately, make falsifiable predictions. Stronger and stronger constraints 
could be placed on the free parameters without ever ruling out the effects, unless, 
of course, there is an indisputable observation of Lorentz invariance violation. It is 
hoped that future developments in causal set theory will allow 'natural' values for the 
phenomenological parameters to be derived. Observations that conclusively rule out 
the models would then be possible. In the meantime the models help us determine 
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where to look for the signature of discreteness. The questions raised by the investiga- 
tions into phenomenology will hopefully also contribute to progress in the underlying 
theory. 

Chapter [2] introduces causal set theory and provides the basic definitions that will 
be required later in the thesis. Chapter [3] discusses the effect of spacetime discreteness 
on the propagation of massive particles. This chapter begins by reviewing the swerves 
diffusion model proposed by Dowker et al. in 131, while the remainder of the chapter 
consists of original work. In Chapter |4]massless particles are considered and found to 
also experience diffusion. Chapter [5] develops the work on massless particles further 
by considering the effect of discreteness on photon polarisation. 

Portions of Chapter |3 together with a large part of the content of Chapter [4] appear 

in 

F. Dowker, L. Philpott, and R. Sorkin. Energy-momentum diffusion from 
spacetime discreteness. Phys. Rev., D79T24047, 2009, 0810.5591. 

Section l331 appears in 

L. Philpott. Particle simulations in causal set theory. Class. Quantum Grav., 
27:042001, 2010, 0911.5595. 

The majority of the work in Chapter [5] appears in 

C. Contaldi, F. Dowker, and L. Philpott. Polarisation Diffusion from Space- 
time Uncertainty, 1001.4545. 

Throughout this thesis Planck units defined by c = h = G = 1 will be used. Unless 
otherwise specified Boltzmann's constant is also set to one: fee = 1. 



Chapter 2 



A BRIEF INTRODUCTION TO CAUSAL 
SET THEORY 



... in a discrete manifold the principle of metric relations is already con- 
tained in the concept of the manifold, but in a continuous one it must 
come from something else. Therefore, either the reality underlying space 
must form a discrete manifold, or the basis for the metric relations must be 
sought outside it, in binding forces acting upon it. 

Riemann (translated in Q) 

The continuity of space apparently rests upon sheer assumption unsup- 
ported by any a priori or experimental grounds. 

Whitehead El 

Given that both general relativity and quantum theory are based on continuum 
spacetimes, it may seem slightly strange that spacetime discreteness is invoked in at- 
tempts to reconcile the two. Why throw away the one thing these theories have in 
common? In brief, the singularities and associated infinite curvatures of general rel- 
ativity could indicate that the continuum is not a good description of spacetime on 
very small scales; in quantum field theory discreteness would give a natural cutoff 
and allow issues with renormalisation to be avoided. 

Discrete spacetime is by no means a new idea. Riemann (as quoted above) points 
out the elegance of a discrete manifold: a discrete manifold actually contains more 
information than a continuum. Einstein, in view of the discrete nature of matter, sug- 
gested that "the continuum of the present theory contains too great a manifold of pos- 
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sibilities".[J Schrodinger also discusses the subject, stating that the continuum "may 
very well turn out to be out of place for physical space and physical time" 0. A 
comprehensive review of the history of discrete spacetime would fill a thesis in and of 
itself; for an interesting discussion of early ideas see [8J. One further historical aspect 
should be mentioned briefly: the first Lorentz invariant 'quantised spacetime' was de- 
scribed by Snyder in 1946 0. It is, however, a rather different approach to the one 
discussed here. 

A considerable difficulty in formulating a discrete theory lies in the fact that many 
of the mathematical tools we are accustomed to using in the description of physics 
rely on the existence of the continuum. What is the analogue of a differential equation 
of motion in a discrete spacetime? Developing a discrete theory of spacetime involves 
not only new physics, but a new language to describe that physics. 

Causal set theory, first proposed by Bombelli et al. in 1987 [10], is a discrete, Lorentz 
invariant approach to quantum gravity. Although there are a number of approaches 
to quantum gravity that are in some sense discrete (Regge calculus, causal dynamical 
triangulations, loop quantum gravity, to name a few), the exact nature of the discrete- 
ness varies greatly. Causal set theory draws its primary motivation from theorems 
proved by Malament ]11| and Levichev |12|, and based on the results of Hawking, 
King, and McCarthy [I3|: for a past and future distinguishing spacetime, the causal 
ordering together with a four-dimensional volume element provides sufficient infor- 
mation to determine all metric, topological, and differentiable structure of a spacetime. 
If spacetime is discrete, this result translates to (in a phrase coined by Sorkin): 

order + number = geometry. 

Number here is just the equivalent of volume - in a discrete spacetime the volume of 
a region is determined by counting the elements within it. In other words, a collection 
of elements endowed with a causal ordering should be all that is needed to describe 
all the complexities of spacetime to the discreteness scale. Early proposals for a theory 
based on discrete spacetime, endowed only with a causal ordering of events, were 
made independently by both Myrheim [14| and 't Hooft 1115] in 1978. These theories 
remained apparently undeveloped until causal set theory was put forward in 1987. 
For comprehensive reviews of the field see, for example, |[T6l[T7l[T8l[T9ll20l . 

1 Einstein to Walter Dallenbach, November 1916, excerpts published in |6|. 
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2.1 What is a causal set? 

Specifically, a causal set is a set C endowed with a binary relation -< satisfying: 

1. transitivity: if x -< y and y < z then x < z, Vx, y,z G C; 

2. reflexivity: i ^ i, Vi 6 C; 

3. acyclicity: if a; -< y and y -< x then x = y, Vx, y G C; 

4. local finiteness: Vx, z G C the set {y | x -< y -< z} of elements is finite. 

The relation -< is commonly called 'precedes'. This definition seems a little more in- 
tuitive if one considers, for a moment, a standard continuum Lorentzian manifold, 
M. The usual causal future relation, J + , satisfies the first three of the above points. 
The causal future J + (x), x G M., is the union of x with the set of all y G M. that can 
be reached from x by a future-directed, non-spacelike curve in A4. This relation is 
clearly transitive, reflexive (x G J + (x)) and acyclic (provided the spacetime has no 
closed causal curves). The condition that the causal future relation does not satisfy is 
local finiteness. Local finiteness provides the discreteness of the theory - in any causal 
interval of the causal set there are only finitely many elements. 

Causal sets can be visualised using Hasse diagrams (see Figure l2T]) . Here time runs 
vertically, elements are drawn as points and relations as lines. Relations implied by 
transitivity are omitted to avoid over-complicating the diagram. Even so, visualising a 
causal set using a Hasse diagram is only feasible for very small causal sets. Although 
Hasse diagrams will not be used in the remainder of this thesis, Figure 12.11 provides 
useful examples for the following important definitions. 

Let C be a causal set. 

1. A chain is a totally ordered subset of C, e.g. a -< /3 -< ( in Figure l2Tl 

2. A longest chain between two elements x, y G C is a chain whose length is longest 
amongst chains between those endpoints, e.g. a ^ (3 ^ 5 ^ ( is a longest chain 
between a and £. There may be more than one longest chain between two ele- 
ments. The length (i.e. number of steps) of the longest chain between elements 
x, y G C will be denoted d(x, y). For example d(a, () = 3. 
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time 




Figure 2.1: An example Hasse diagram of a 7 element causal set. 

3. A link is an irreducible relation: elements x and y are linked if and only if 
d(x, y) = 1, e.g. a -< f3 is a link. If two elements x, y £ C are linked, it will 
be denoted x -<* y. 

4. A path is a chain consisting of links, e.g. a -<* f3 -<* 7. 

2.2 The causal set hypothesis 

Above I provided a mathematical definition of a causal set. Causal set theory is based 
on the hypothesis that spacetime is a causal set (or, indeed, a quantum sum-over- 
causal sets). The observed continuum Lorentzian manifold, it is assumed, arises as an 
approximation to an underlying causal set. The partial order gives rise to the causal 
ordering of events in the approximating continuum spacetime, and the number of ele- 
ments comprising a spacetime region gives the volume of that region in fundamental 
units. The fundamental volume unit is expected to be of the order of the Planck vol- 
ume (where necessary in the later chapters of this thesis, it will be assumed that the 
fundamental volume is equal to the Planck volume). Note that any continuum man- 
ifold that approximates a causal set is necessarily Lorentzian, as only a Lorentzian 
metric can give rise to a partial order on spacetime points. 

The central, and as yet unproved, conjecture of causal set theory is that if two 
continuum manifolds are in some sense a 'good approximation' to a given causal set, 
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Figure 2.2: A diamond lattice no longer appears uniform when viewed in a boosted frame. 



then those two continuum manifolds should themselves be approximately the same. 
The idea here is simply that although multiple manifolds could approximate a given 
causal set, the manifolds should capture all the significant information from the causal 
set and only vary between themselves on very small scales. It is far from simple to 
mathematically define what is meant by 'approximate' in this situation. Some intuitive 
understanding can be gained by working backwards: constructing causal sets from 
manifolds via a process called 'sprinkling'. 

2.3 Sprinklings 

A causal set can be constructed by selecting points from a Lorentzian manifold. The 
causal relations between the points on the manifold induce the order relations between 
the elements in the causal set. The most obvious way to select points from a manifold 
is to construct a regular lattice, but this does not provide the features needed in causal 
set theory. Consider choosing a frame and selecting points from two dimensional 
Minkowski spacetime in a diamond lattice. If a Lorentz boost is applied the points are 
no longer a good sampling of the manifold, as is shown in Figure 12.21 In the boosted 
frame the number of elements in any region is clearly not a good approximation of the 
volume of that region. 

Consider, instead, selecting points at random from a manifold via a Poisson pro- 
cess in which the probability measure is equal to the spacetime volume measure in 
some fundamental units - this process is called 'sprinkling'. The number of points 
chosen from any region of the manifold will be approximately equal (up to Poisson 
fluctuations) to the volume of that region in fundamental units. If a Lorentz boost is 
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Figure 2.3: A causal set produced by sprinkling into a region of two dimensional Minkowski 
spacetime remains uniform density under a Lorentz boost. The empty regions in the boosted 
case are edge effects due to sprinkling into a finite region. 

applied to points selected this way, the distribution remains of uniform density (see 
Figure l23l) . For a proof that the sprinkling process is Lorentz invariant, see \21\. Sprin- 
klings prove very useful in analytic and numeric studies of causal sets - Chapter[3]will 
make use of them in both ways. 

Earlier in this section, the longest chain between two elements x, y in a causal set 
was defined. For a causal set generated by sprinkling into Minkowski spacetime, the 
length of the longest chain, d(x, y), is found to be proportional to the proper time 
between x and y in the limit of large distances |22|. The constant of proportionality 
is dependent on the spacetime dimension. A longest chain in a causal set is thus 
a close approximation to a timelike geodesic. For a numerical investigation of this 
correspondence in conformally flat spacetimes see 11231 . 

2.4 Causal set phenomenology 

As discussed in Chapter [T] the majority of current research in quantum gravity phe- 
nomenology focuses on violations of Lorentz invariance. Causal set theory, on the 
other hand, provides a way of investigating Lorentz invariant quantum gravity phe- 
nomenology. The difficulty of course, is that causal set theory is incomplete, and ob- 
taining testable predictions from an incomplete theory is far from simple. Although 
much will have to wait until the development of a quantum dynamics for causal 
sets, the concreteness and simplicity of the causal set hypothesis allow considerable 
progress to be made even at our early state of research. 
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Before moving on in Chapter[3]to the particle phenomenology that forms the basis 
of this thesis, another significant result in causal set phenomenology should be men- 
tioned. As early as 1997 Sorkin [24J suggested that causal set theory could give rise 
to a cosmological constant of the order 1 CP 120 - a number that matches current ob- 
servations (see also 1131261). This prediction for, in fact, a non-constant cosmological 
constant relies on the simple idea that a discrete spacetime leads to fluctuations in vol- 
ume. As discussed above, the sprinkling procedure places N points in a continuum 
volume V with fluctuations of the order y/~N. Although there are questions that re- 
main to be addressed fl27l |28| this model illustrates the powerful results that can be 
obtained from causal set theory. 



Chapter 3 



Swerves 



quare etiam atque etiam paulum inclinare necessest 
corpora; nec plus quam minimum, ne fingere motus 
obliquos videamur et id res vera refutet. 

Lucretius, De Natura Rerum, 11. 243-245 

Wherefore it is necessary that bodies swerve a little perpetually, 

Not more than a little, lest we seem to imagine sideway motions, and the 

facts confute our conjecture O 

General relativity predicts that a free, massive particle will travel along a timelike 
geodesic - a prediction well tested by observations. If the underlying structure of 
spacetime is discrete rather than continuous, particles will no longer be able to travel 
along perfect geodesies. Discreteness, one expects, will introduce some random fluc- 
tuations into a particle's trajectory. Of course, given how well tested general relativity 
is, if such fluctuations exist they must be so small as to have escaped notice. Space- 
time discreteness is likely to occur on the Planck scale, so it is certainly reasonable to 
suppose that an effect on particle trajectories exists that has not yet been observed. 

In |3| Dowker et al. consider this problem in the context of causal set theory. Moti- 
vated by a simple model of particle motion on a causal set (see Section 133) they sug- 
gest that particles travelling through a discrete spacetime will be subject to swerves: 
random fluctuations in momentum. From a simple random walk on R we can derive 
the standard diffusion equation in the continuum limit: similarly, a basic stochastic 
process that gives random fluctuations in momentum is found in the continuum limit 
to lead to a diffusion equation. 

translation courtesy of Dr. R. M. Pollard, who notes that inclinare would be better translated as 'move 
back and forth', which is somewhat less succinct yet even more fitting for our context. 
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The swerves diffusion equation given in [3] does not depend on any specific un- 
derlying particle model - nor, in fact, does it explicitly depend on the hypothesis that 
spacetime is a causal set. The power of the swerves diffusion equation lies in it be- 
ing the unique Markovian Poincare invariant diffusion equation on the state space. 
Any underlying process that gives small random fluctuations in momentum will be 
described by this equation in the continuum limit. 

A full derivation of the swerves diffusion equation, first stated in 0, is given in 
Section I3T21 The diffusion equation, as derived by Dowker et al. in terms of the proper 
time of the particle, is not the most useful form. To make contact with observations an 
equation in terms of an observable cosmic time is needed. In [3] Dowker et al. suggest 
the form of such an equation in the special case of a homogeneous distribution of 
particles. Here this work is extended and a full derivation of the inhomogeneous case 
is given in Section 1331 

I begin, however, by investigating models for particle propagation on causal sets 
more fully. Three models, including the model first described in O, are discussed 
in Section 13.11 below. The connection between fundamental models and the diffu- 
sion equation is investigated in Section l3~4l where it is shown that the swerves model 
on the causal set gives rise to a finite nonzero diffusion constant. The work in Sec- 
tions |3T]3i2l33] appears in [29| and together with Section l3~4l is work done in collabo- 
ration with Fay Dowker and Rafael Sorkin. Section 1531 discusses the observable con- 
sequences of the swerves diffusion equation and reviews the bounds placed on the 
diffusion constant in ©1301. Lastly, the models introduced in Section I37T1 are investi- 
gated numerically in Section 13.61 and the behaviour is compared with that predicted 
by the swerves equation. The work in Section 1331 appears in II3TH . 
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3.1 Models for massive particle propagation on a causal set 

What is a particle in causal set theory? It must be admitted that this is not a question 
that can currently be answered. In general relativity matter determines spacetime cur- 
vature. In the context of a causal set this curvature must be encoded in the structure of 
the causal links - perhaps a particle is some pattern or knot of causal links? Of course, 
if the mass of the particle is a manifestation of a tangle of links between causal set 
elements, what about the other particle properties, the charge or spin? Approaching 
from the direction of quantum field theory we could assign amplitudes for a particle 
to be 'on' a particular causal set element and to move from one element to another. 
For work on developing quantum field theory on causal sets, see Il32ll33| . 

An understanding of particles in causal set theory will really require a theory of 
quantum causal set dynamics (QCSD). If QCSD takes the form of a 'sum-over-causal 
sets', will a particle trajectory be some sum over all trajectories in all possible causal 
sets? Even without a theory of QCSD an intuition into the effect of discreteness can 
be gained by considering simple models. Massive particles will be considered as clas- 
sical point particles located 'on' a causal set element, a particle trajectory will con- 
sist of jumps from one element to another. The models discussed in this section are 
clearly not physically realistic, but they do capture the important aspects of the causal 
set approach and therefore provide one very useful way of investigating observable 
consequences of discrete spacetime. The models introduced in this section will be 
investigated numerically in Section 1331 

3.1.1 The swerves model 

This model was introduced in J3]| . 

Model 1 

Construct a causal set C by sprinkling into Minkowski spacetime. A massive particle 
trajectory is taken to be a chain of elements e n in the causal set, i.e. a linearly ordered 
subset of C. It is assumed that the trajectory's past determines its future, but that only 
a certain proper time tj into the past is relevant. When this 'forgetting time' Tt is small 
the process should be approximately Markovian. t/ is also assumed to be much larger 
than the discreteness scale. The particle trajectory is constructed iteratively. Suppose 
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the particle is currently located 'on' an element e n , with a four-momentum p n . The 
next element, e n+ \ is chosen such that 

• e n+ i is in the causal future of e n and within a proper time tj of e„, 

• the momentum change \p n +i — Pn\ is minimised. 

Here the momentum p n +i is defined to be proportional to the vector between e n and 
e n+ i, and thus this model relies on knowledge of the embedding of the causal set into 
Minkowski spacetime. The constant of proportionality is chosen such that the mass of 
the particle remains fixed, i.e p 2 = —m 2 . This method for determining the trajectory is 
illustrated in Figure 13.11 where the frame has been chosen to be the rest-frame of the 
particle at e n . Heuristically, the requirement that the momentum change be minimised 
creates a trajectory that stays as straight as possible at each step. If we work in the rest- 
frame of the particle, this method chooses the next element to be as close as possible to 
the t axis within the forgetting time. At each step there is a small random fluctuation 
in the momentum of the particle. There is no direct dependence on the mass of the 
particle in this model. Any nonzero mass will result in the same trajectory, although 
we could imagine that the forgetting time t/ is mass dependent. 
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This model is clearly not a realistic fundamental law of motion for a particle on a 
causal set. It requires information not contained in the causal set - the embedding of 
the elements into the approximating Minkowski spacetime. This issue is overcome in 
the 'intrinsic' models described below. All these models treat particles as classical and 
zero size. 

3.1.2 Intrinsic models 

For a model of particle propagation to be intrinsic to a causal set, it cannot refer to 
a continuum forgetting time Tt. Instead, a 'forgetting number' (an integer nt) must 
be defined. The forgetting number can be interpreted as the number of discrete steps 
into the past of the trajectory that are relevant in determining the future trajectory 
(although in Model 2 below this number is in fact 2nA If the embedding into a con- 
tinuum spacetime was known, nj could be roughly written nj = Tf/d p i, where d p i is 
the discreteness scale. It is important to note, however, that the models are defined on 
a general causal set and it is not assumed here that the underlying causal set can be 
faithfully embedded into a 'reasonable' spacetime. 

In defining these intrinsic models there is no reference to the mass of the particle. 
As above, the forgetting number could be taken to be mass dependent. The models 
discussed cannot, however, be considered as models of particle propagation for mass- 
less particles - these models give trajectories whose long time behaviour approximates 
a timelike geodesic, rather than the null geodesic appropriate to massless particles. 

Two slightly different intrinsic models will be described in this section, to give an 
idea of the wealth of possibilities available. First recall some causal set definitions 
from Chapter [2) 

• d(a, b) denotes the length of the longest chain between two elements a and b, 

• a path is a chain consisting of links a -<* b, i.e. d(a, b) = 1. 

Model 2 

A massive particle trajectory is taken to be a chain of elements . . . e n _2 -< e n _i -< e n 

Given a partial particle trajectory . . . e n _i, e n the next element e n+ \ is chosen such that 

• d(e n -i,e n+ i) < 2n { , 
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• d(e n , e n+ i) is maximised subject to d(e n , e n +i) <Tif, 
(see Figure 3.2(a)| ). These requirements do not guarantee the existence of a unique 



e n+ \. There will, however, almost surely be finitely many eligible elements and the 
trajectory can be constructed by choosing an element uniformly at random from these. 
Note that this model is slightly different from the first intrinsic model presented in [29J 
(where equalities in the above conditions were given). Model 1 of B29II does not guar- 
antee the existence of an e n+ \ under reasonable conditions. 

Under this model the particle trajectory should swerve a little, but remain approx- 
imately straight so long as nj is large. It is easiest to see this if we consider the 'ideal' 
case where there exist e n _i, e n , e n+ \ such that d(e n -i,e n ) = n/ and d(e n , e n+ \) = nf. 
The elements have been chosen such that d(e n -\, e n +i) < 2n/; since d(e n _i, e n ) = rif 
and d(e n , e n+ i ) = nf we know there exists some chain between e n _ i and e n+ i of length 
2rif and hence the equality d(e n _i, e n+ i) = 2nj must hold. In other words, the chain 
we have chosen between e n _i and e n+ \ is a longest chain, and thus the trajectory is 
approximately geodesic over {e n _i, e n _|_i} segments. In practice it will not always be 
possible to choose an e n+ \ such that d(e n , e n+ \) = rif and d(e n -i, e n +i) = 2nj, neces- 
sitating the inequalities. 

In this model the trajectory can be considered as composed of just the elements 
. . . e n _i, e n , e n+ i ... or of the 'filled in chain' consisting of a (randomly chosen) longest 
chain (of length < nf) between e n _i and e n , another between e n and e n+ \, and so 
on. Variations on this model can also be constructed. One possibility is to choose the 
forgetting number at random at each step from a distribution with a mean rif and 
some fixed variance. 

Model 3 

The trajectory is constructed as a path in this model, i.e. d(e n ,e n+ i) = 1 for any 
e n , e n+ \. Given a partial particle trajectory . . . e n _ n/ , . . . ,e n _i,e n the next element e n+ \ 
is chosen such that 

• d(e n ,e n+ i) = 1, 

• d(e n - nf , e n+ i) + . . . + d(e n -i,e n+ i) + d(e n , e n+ i) is minimised, 



(see Figure 3.2(b)| ). Note that this minimisation does not necessarily yield a unique 



e n+ i, in which case we construct the trajectory by choosing an element uniformly at 
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d(n - l,n + l) < 2n f 



d(n + l,n + 2) < n/ 




tZ(n, n + 1) < ra/ 



d(n - 1 , n) < 71/ 



e. 



d(n,n + 2) < 2n } 



d{n-l,n + 2) > 3 




d(n,n + l) = 1 



d(n — l,n) = 1 



d(n + l,n + 2) = 1 



(a) Model 2. 



(b) Model 3. 



Figure 3.2: Trajectories constructed using intrinsic models. 



random from those eligible. Also, if the trajectory has length less than rij the minimi- 



sation is done over all elements available. 

In this model each element is linked to the previous, i.e. d(e n _i, e n ) = 1, and thus 
we know there exists a chain (our trajectory) of length n$ between e n _ n/ and e n . The 
longest chain length, d(e n - nf , e n ), must therefore be greater than or equal to n/. If we 
choose e n to minimise d(e n - nf , e n ) we ask that the trajectory be as close as possible to 
geodesic between e n _ n/ and e n while fulfilling d(e n -i, e n ) = 1. Minimizing the sum 
of the partial lengths distributes the geodesic property along the path. 
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3.2 The swerves diffusion equation 

Although the intrinsic models described in Section 13.11 are defined on any causal set, 
the process is really only of interest when the causal set is a good approximation to 
Minkowski spacetime (or another physically relevant spacetime, such as Friedmann- 
Robertson- Walker). The universe we live in is approximately Minkowskian, and thus 
even though the behaviour of particles on nonembeddable causal sets may be very in- 
teresting, it is of little relevance when investigating the observable phenomenology of 
discrete spacetime. If it is assumed that the causal set can be produced by sprinkling 
into Minkowski spacetime, the particle models can be thought of as defining piece- 
wise linear curves in Minkowski spacetime. The particle's momentum - which cannot 
(currently) be defined in a manner intrinsic to the causal set - is then defined every- 
where on the curve except at the vertices. It is clear here that in the continuum limit 
the intrinsic models lead to random fluctuations of momentum (defined, as usual, in 
Minkowski spacetime). 

Rather than choosing a specific particle model on a causal set and investigating its 
continuum behaviour, an entirely general situation will be considered. Motivated by 
the models for particle trajectories on a causal set, it is assumed that there exists some 
underlying Markovian, Poincare invariant process that causes random fluctuations in 
momentum (and consequently position). The process takes place in the proper time 
of the particle. Note that the models described above are approximately Markovian 
when the forgetting time (number) is small. If the discreteness scale is of the order of 
the Planck scale, the forgetting time can be many orders of magnitude greater than the 
discreteness scale and yet still small compared to the trajectory length. 

3.2.1 Stochastic evolution on a manifold of states 

The swerves equation is obtained using the general formalism of |34| for stochastic 
evolution on a manifold of states. Chapters H] and [5] also draw on this framework, and 
thus for completeness the main points in this work will be discussed here. 

Consider a 'mesoscopic' process occurring on a manifold of states M.. The current 
state of the system is described by a probability density p. It is assumed that the future 
of the physical system depends only on its present state and not its past, i.e. the evo- 
lution is a Markov process. It is also assumed that the system traces a continuous path 
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through M . Such a process can be described by a linear, first order in time equation 



where p = p[x, T), the probability density for the system to be in a state x G M at time 
T, is a scalar density in M. The indices A,B = 1,..., dim(M). K AB = K BA , B A , C are 
functions of x (and possibly T). For extensive justification of this assumption see |34|. 
Note that if Equation l3. 1 1 were to contain higher order terms it would allow unphysical 
negative probability densities, p < 0, proof of this is given in Appendix A of [34 J. 

The requirement that the probability density p be nonnegative leads to a constraint 
on K AB in the above, general equation. Suppose the initial distribution, in some 
neighbourhood of the origin, is given by 



where Z must be positive by the requirement p > 0. Clearly p(0, 0) = and in order for 
p(0, T) > we must have p(0, 0) > 0. From Equation [3j] we can conclude Z AB K AB > 
0. Suppose Zab = zazb- This implies zazbK ab > 0, i.e. K AB is a (symmetric) 
positive semidefinite matrix (a matrix A is positive semidefinite if x T Ax > Vx G M n ). 
We can also impose the conservation of probability 



With p given by Equation 13.11 integration by parts implies C = 8a9bK ab — d A B A . 
Substituting this expression back into Equation 13.11 gives 



J- = K AB d A d B p + B A d AP - Cp , 



(3.1) 




(3.2) 




(3.3) 



dp 
dT 



K AB d A d BP + B A d AP - {d A d B K AB - 8 A B A ) p 

d A d B (K AB p) + d A {B A P ) - 2d A d B K AB P - 2d B K AB d AP 

d A d B (K AB p) + d A (B A p) - 2d A (d B K AB p). 



(3.4) 



Equation 13.11 can therefore be written in terms of a current, J A , and a continuity 
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equation: 



J A = -d B (K AB p)+v A p, (3.5) 
-d A J A , (3.6) 



9 P :1 r.i 



dT 

where u A = 2d B K AB - B A . Note that J A is a vector density (see pg. 124 of (33) and 
the vector density transformation implies that while K AB is a tensor, v A is not. The 
above equations can be expressed in terms of a true vector u A if the density of states, 
n, is introduced. The vector u A is defined by 

u A = v A - d B K AB - K AB d B s , (3.7) 

where s is the entropy scalar 

s = k B log n . (3.8) 

Here k B is the Boltzmann constant (note that in the sections that follow k B = 1). 
Equations I3.5l and l3.6l can then be written in the form 

^ = d A (K AB D B p - u A p) , (3.9) 

where 

Dap = d A p - pd A s . (3.10) 

In the sections that follow, it will often be useful to use the following relations for K AB 

and v A : 

K AB = to ( A t A f) , (3.1D 
AT^o\ 2AT / ' v ' 



A 



v A = lim ( —— ) . (3.12) 

AT^O \ AT / v ' 

Here < > denotes the expectation value. These equations give a way of relating the 
abstract objects K AB and v A to the basic properties of the physical stochastic process 
- the 'spatial' step Ax and the time step AT. Equations 13.111 and 13.121 are derived in 
Appendix B of |34| . 



§3.2 The swerves diffusion equation 



32 



3.2.2 Diffusion in proper time 

The state space, M., of a swerving particle of mass m is M. = M 4 x H 3 , where H 3 is 
the mass shell. A point on M. thus represents a position in M 4 and a momentum in 
H 3 . The coordinates on M 4 are the usual Cartesians {x^ 1 }, \i = 0,1,2,3 and indices 
are raised and lowered with rj^, the Minkowski metric. The spatial coordinates on 
M 4 will be written as {x 1 }. Cartesian coordinates in momentum space are p^, where 
Pn is subject to the constraint that it lies on the mass shell, i.e. the hyperboloid H 3 in 
momentum space defined by p^+m 2 = 0. p° = E is the energy (taken to be positive) 
and p = \Jp\ + p\ + p\ is the norm of the three momentum. The three coordinates on 
H 3 will be written abstractly as p a . Coordinates on M. are denoted collectively as 
X A = {x^,p a } and in what follows capital letters A, B will be used to indicate general 
indices on M; fj,, v = 0, 1, 2, 3 are indices on M 4 ; i, j = 1,2,3 are spatial indices on M 4 ; 
a, b = 1, 2, 3 are indices on H 3 . 

The metric on M is the product of the Minkowski metric r]^ on M 4 and the 
Lobachevski metric g a b on H 3 . This is the unique Poincare invariant metric on Ai 
(up to an overall constant). The 'density of states', n, plays a role in the formalism 
of |34'|, and by symmetry, it must be proportional to the volume measure on M, so 
n oc where g = det(g ab ). 

As described above, a process that undergoes stochastic evolution on a manifold 
of states, M, in time parameter T, can be described by a current, J A and a continuity 
equation II34H : 



J A 
dT 




d A J A , 



(3.13) 



(3.14) 



or alternatively by the equation 




(3.15) 



where 



u 



A 



= V 



,A 



d B K 



AB 




(3.16) 



To find the diffusion equation for the swerves particle process, it is necessary to deter- 
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mine K AB and u A . The requirement that the equation be Poincare invariant is a very 
stringent condition and allows K AB and u A to be determined up to the choice of one 
constant parameter. 

Consider the process in terms of r, the proper time along the worldline of the 
particle. From Equation l3.11l in Section l3.2.1l 

„ AH , /AX A AX B \ 

K AB = lim ( ) . (3.17) 



At^-o+ \ 2Ar 

Recall K AB is a symmetric, positive semidefinite matrix, and here is also required to 
be Poincare invariant. Looking first at the spacetime component of this matrix, 

K~- Mm (^1). (3 . 18) 
Ar^o+\ 2Ar / v ' 

Ax 11 = ip^Ar at every step of the process and so K^ v = ^\\m/^ T ^Qp^p v At = 0. 
Given K^ u = 0, K^ a = is required by the condition that K AB be positive semidef- 
inite. To see this, first recall that a matrix is positive semi-definite if all the principal 
minors of the matrix are nonnegative. Consider a symmetric 2x2 matrix of the form 



(3.19) 



To be positive semi-definite the matrix must have —a? > 0, i.e. a = 0. Similarly, 
beginning with 

/ K^ a 




K AB =1 I , (3.20) 



we can consider all the 2x2 principal minors and determine K^ a = K afl = 0. 

The only Lorentz invariant tensor on H 3 is proportional to the metric, g ab , and the 
coefficient is independent of x^ 1 by translation invariance, thus K ab = kg ab , where 
k > is a constant. This gives 



K AB = I ° ° I . (3.21) 



kg 



ab 
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To determine u A first recall Equation l3.12l 

v a = lim / ^1) . ( 3. 22) 
Ar^o+ \ At I 

The spacetime component 

/ Ar^\ 

if = lim ( — — ) , (3.23) 
Ar^o+ \ At / 

is simply = p^/m. The components of the true vector u M are equal to since 
Rf iA = _ There 

is no Lorentz invariant vector on H 3 and so u a = 0, giving 

u A = (p^/m,0) . (3.24) 



The proper time diffusion equation can now be written down from Equation l3.15l 

£-*fl.(/V»(^))->%*. (3.25, 

This can be seen to be equivalent to Equation 1 of 131 if a scalar p = p T /^/g is 
defined: 

^- = kV 2 H p--p^d tl p, (3.26) 
or m 



where V 2 H is the Laplacian on H 3 . 
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3.3 Diffusion in cosmic time 

Given an initial distribution of particles, for instance from an astronomical source, the 
equation derived above is not very useful for predicting the results of observations. 
Even if particles all leave the source at the same time with the same momentum, the 
momentum variation induced by the swerves will result in particles arriving after dif- 
ferent proper times and at different observatory times. The proper time that elapses 
along the particles' worldlines from source to detector is not observable. To compare 
the swerves model with experiment and observation it is necessary to describe the 
evolution of the distribution in time in the rest frame of our detector, which time will 
be referred to as cosmic time. One may ask why the equation was derived in terms 
of proper time if it is more useful to obtain a cosmic time equation. The reason be- 
comes clear if one considers that when working with proper time, the explicit Lorentz 
invariance makes it simple to determine K AB and u A . The result can then be easily 
transformed to a specific frame, as shown here. 

A first step in this direction was to look at the nonrelativistic limit of the proper 
time diffusion equation, when proper time and cosmic time are comparable. The non- 
relativistic limit in fact proves sufficient to place very strong bounds on the value of 
the diffusion constant and severely limit any observable effects (see |3| and |30| and 
Section|33]). 

In the fully relativistic case, Dowker et al. wrote down the diffusion equation in 
terms of cosmic time for the special case of an initially spatially homogeneous distri- 
bution H3j. The derivation of the cosmic time evolution equation for the general case 
of a spatially inhomogeneous distribution will be given here. 

The conversion between proper time and cosmic time is possible because both are 
good time parameters along all possible particle worldlines, which are causal. If the 
diffusion process is visualised as a collection of worldlines through spacetime and 
momentum space, both cosmic time, t = x°, in our chosen frame and proper time r 
increase monotonically along each trajectory. Assume that the particle starts at param- 
eter r = and cosmic time t = 0. Proper time can be added to the state space and the 
process is represented by flowlines in M! = M 4 xtfxl (see Figure l33l) . Along each 
flowline, both r and t are good time parameters. The proper time diffusion equation, 
Equation l3.25[ describes the evolution of the distribution on constant r hypersurfaces 
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Figure 3.3: Particle trajectories as flowlines in Ad' = M 4 x H 3 x K (where the momentum and 
two spatial dimensions have been suppressed). 

in M! . What is needed is a diffusion equation for evolution of the distribution on 
constant t hypersurfaces integrated over all proper times. 

First, working in the larger space M! , define a new current component 

J T (t,x\p a ,r) = p T . (3.27) 

If coordinates on this extended space M! = A4 x R, are denoted by X a = {X a ,t} 
then the continuity equation (|3.6[) can be written 

d a J a = . (3.28) 

Using Equation l|3.5[) (and still treating r as the time parameter) the t component of the 
current can be expressed in terms of J T (equivalently p T ). 

j\t,x\p a ,T) = -d B {K tB J T ) +v t J T 

= v l r 

= lJ\ (3-29) 



where 7 = dt/dr is the usual relativistic gamma factor. 
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The remaining components of the current can now be written in terms of J*. The 
spatial components are: 



J i (t,x i ,p a ,r) = -d B (K iB J T ) + v* J T 

= v\r 
p 1 j 



i jt 

(3.30) 

m 7 

In the case of the p components the algebra is simpler if Equation 13.51 is expressed in 
the form (see Equation l3.15|) 

jA = _ R AB n d JP\ + pu A^ (3 31) 



n 



and so 



J a (t,x i ,p a ,T) = -kg ab nd b ' Jl 



= -k^y/gdb^-^j • (3.32) 

The metric g ab that appears here is the Lobachevski metric on H 3 . 

Since r is unobservable we need to integrate J over r: the current J describes 
the flow through a region AtAxAp at any point (t,x l ,p a ,r); what is needed is the 
cumulative flow through AtAxAp at any given (t,x l ,p a ) over all proper times. The 
integrated current will be denoted J. Integrating the t component of the current over 
proper time from zero to infinity gives the probability density on a hypersurface of 
constant t: 

p t = J\x\p a ,t) 

J 1 dr. (3.33) 



33.3 Diffusion in cosmic time 38 



The components of the new current can be written: 



J l (x\ P a ,t) = I J l dT 



dr 



7717 
m 7 
m 7 ' 



J a (x\p a ,t) = j J a dr 
= -kg ab nd b 



(3.34) 



= -kg ab nd b (^-\ . (335) 
\7>V 

Integrating the continuity equation over r gives 

[JX + dtJ* + diJ 1 + 8 a J a = . (3.36) 

J T | T=0 is zero for all t > and J r tends to zero as r goes to infinity for finite t. So for 
all t > 

d t J t + d i J l + d a J a = Q. (3.37) 

Finally substituting Equations 13341 and l335l into Equation [337] and recalling J* = p t , 
gives the cosmic-time diffusion equation 

t = -£*» + **(^(^))- <338) 



This is a powerful phenomenological model because it depends on only one parame- 
ter, the diffusion constant k. Data can therefore strongly constrain k. 
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3.4 Finiteness of the diffusion constant 

The swerves diffusion equation derived in Section 13.21 is independent of the details 
of the underlying particle model. To make clearer how the underlying model on the 
causal set gives rise to a diffusion equation on M 4 x H 3 the swerves model (Model 1) 
will be investigated more fully. Specifically in this section it will be shown that Model 
1 gives rise to a finite nonzero diffusion constant in the macroscopic limit. 

Recall that a process undergoing stochastic evolution on a manifold of states can 
be described by 

J A = -d B {K AB p) +v A p, (3.39) 
% = ~d A J A , (3.40) 

where K AB is a symmetric, positive semi-definite matrix given by 

K". Urn ( AX *i XB ). (3.41, 

In Section I3T21 it was shown that due to the constraint of Lorentz invariance K AB has 
the form: 

/ 

kg ab 

where k is some constant, irrespective of the underlying model. It is not immediately 
clear, however, that the constant of proportionality, k, need be finite and nonzero for 
any particular model. For the specific case of the swerves model, the component K ab 
is indeed nonzero and finite (and thus K AB is given by Equation l3.42|) . Demonstrating 
this is the task of this section. 

In the swerves model there are two length scales: the discreteness scale d v \, and the 
forgetting time Tf, where t/ >> d p i. The trajectory length (or total proper time of the 
particle) must also be much greater than Tf. It is only in the continuum limit d p i — > 
that the model can be described by a diffusion equation. 

Suppose, for the swerves model, that K ab = kg ab . From Equation 13.411 we also 
have 

K*= lim (3.43) 



K AB =1 I , (3.42) 



Ar^o+ \ 2Ar 
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Here Ap a is the change in the momentum component p a in a single swerves step and 
At is the change in proper time for that step. Note that in the swerves model At is not 
constant for each step: it is a random variable in the range < At < t/. We wish to 
determine K ab as d p i — > 0. First note the following two lemmas for the swerves model: 

Lemma 3.4.1. Ifd p i -> and r f is fixed then / Ap 2 \ ->■ 0. 

Lemma 3.4.2. Ifd p i is fixed and Tf — > then (^ A 2A^ P °°- 

Section [3.4.11 below contains the proofs of these lemmas. From these two lemmas 
it is clear that as we take the continuum limit d p i — > we can also take rt — > (and 
thus At — > 0) in such a way that / A ^^ p \ remains fixed at some finite nonzero value. 
Thus, if K ab = kg ab the diffusion constant k is finite and nonzero, i.e. the swerves 
model results in the diffusion equation. 

3.4.1 Proofs of Lemmas 13.4.11 and 13.4.21 
Hyperbolic coordinates 

For ease of calculation, hyperbolic coordinates (77, £, 9, 4>) in four dimensional Minkowski 
spacetime can be defined: 

= t 2 -r 2 , (3.44) 

f = tanbT 1 Q , (3.45) 

9 = arccos (-) , (3.46) 



-r 

arctan 



(I) ' (3 ' 47) 



where r = \J x 1 + y 2 + z 2 . These coordinates are illustrated in Figure 13.41 77 defines 
hyperbolae and £ defines radial directions, varying from £ = on the i-axis to £ = 00 
at t = r. 



Calculating volumes in hyperbolic coordinates 

Before addressing the Lemmas I3.4.1l and l3.4.2l it is useful to calculate the volume of the 
shaded region shown in Figure 13.41 the volume of 4-dimensional Minkowski space- 
time such that < 77 < rj and < £ < The first step is to determine the vol- 
ume element in coordinates (rj, £, 9, 4>). It is easiest to begin with spherical coordinates 
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(t,r,6,(j)) defined by 

t = t, (3.48) 

x = r sin 9 cos <p , (3.49) 

y = r sin 9 sin cf) , (3.50) 

2 = rcosO. (3.51) 

The spherical coordinates can then be expressed in terms of hyperbolic coordinates 

t = rj cosh £ , (3.52) 

r = ??sinh£, (3.53) 

= 9, (3.54) 

= 0. (3.55) 
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The Jacobian determinant for this change of variables is 
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The volume element is thus given by 



(3.56) 



dt dx dy dz = r 2 sin 9 dt dr d9 d<j> 

= rf sinh 2 £ sin 9 dr] d£ d9 1 



(3.57) 



The volume of the shaded region is 

r-2n fir /•£' rrf 
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rf sinh £ sin 9 dr] dt; d9 1 
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/•£' rrf 
4tt / rf sinh 2 £ dr] d£ 
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,4 / f ' sinh 2f ' 



TTT] 



2 + 4 



(3.58) 



This result will be used repeatedly in the calculations below. To prove Lemmas 13.4.11 
and !3.4.2l the following property of sprinklings is required: if points are sprinkled with 
a density p into a region, the probability that n points are sprinkled into a volume V is 
given by the Poisson distribution 



Pr, 



(pV) n e 



(3.59) 
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For example, the probability that no points are sprinkled into the shaded region in 
Figure l3~4l is simply 

Pr p = e- pv 



( ,4/ e smh2? \\ 
= expl-pvrr? (~^ + 4 ))■ ( 3 - 60 ) 

To prove Lemmas 13.4.1 1 and 13.4.21 it is necessary to calculate the expectation value 
A 2A^ P )• m ^ act/ srnce ^ i s known that K a6 oc g ab it suffices to calculate a single 
component, 

Working in terms of polar coordinates {p, 0,(j)} on H 3 , the metric g ab is 



2 

ds 2 = m dp 2 + p 2 (dfl 2 + sin 2 eta 2 ) . (3.61) 
m 2 + p l 

Consider a swerves trajectory in the rest frame of the particle after n steps (see Fig- 
ure [37TJ . In this frame the metric component g pp = m = 1 and thus K pp = k. Since 
the three momentum p n is zero, the change in momentum in the next step is simply 
p n +i. Working in Cartesian coordinates for a moment, the momentum at step e n+ \ 
is proportional to the coordinates X^ of the point e n+ \ in this frame: p^+i = ^X^, 
where rj is the proper time between e n and e n+ \. Thus the magnitude of the three 
momentum is given by 



Pn+l = \JPi+P2+pI 



^xf + xl + xl 

m 
— r 

V 



= msinh£. (3.62) 

The expectation value 

where P = P ( ^^ j ^ s * ne probability of the particular value occurring. The proba- 
bility P is the probability that the element e n+ \ has (hyperbolic) coordinates 
rf <rj < rf + drj, £' < £ < £' + d£. Recall that in the swerves model e n+ i is chosen 



33.4 Finiteness of the diffusion constant 44 



such that 77 < Tf and the momentum change is minimised. P is thus the probability 
that there are no points in the region {77 < Tf, £ < £'} and there is a point in the region 
{7/ < 77 < r?' + ^77, £' < £ < £' + The probability that there are no points in the 
region {77 < 17, £ < £'} is simply 



exp I -pttt; I ~ + ==±- ) ) , (3.64) 



f sinh2£' 



2 



from Equation l3.60l above. The probability that there is a point in the region {7/ < 77 < 

77' + dr], £' < £ < g + d£} is 

4^7? ,3 sinh 2 £ , cZ?7(i£, (3.65) 

from the integral in Equation 13.581 and Equation 13.591 (the exponential term can be 
neglected as the volume is very small). The expectation value is therefore (dropping 
the primes on £ and 77) 

/ApAp\ H f°° 2 2 . , 4(t / J £,sinh2A\ 

( 2Ar > = / / 2vrpm 77 sinh £exp I -pvrTy I — - + — - — I J dr]d£ 

= \,pm\) ^ sinh 4 i exp (-prf (- 1 + ^) ) • (3.66) 
Integrating by parts, this reduces to 

<^> - 17 r h (-! + ^)) « ■ <^ 

The proofs of the Lemmas now follow without difficulty. 



Proof of Lemma [3AT1 

If we take the discreteness scale dpi to zero, while keeping the forgetting time, Tf, fixed, then 

Taking d p i — > is equivalent to taking the sprinkling density p to infinity. To deter- 
mine the effect of p — » 00 on Equation l3.67l the integral can be divided into two parts. 
There exists some £ such that, for all £ > |, cosh £ < 2 sinh £. Let 

(^H+, 2 , (3,8, 
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where 
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(3.69) 



Ji can be written 
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(3.70) 



Thus 



lim J x = . 

p— ¥00 



Now consider the integral I2 and the standard result 



(3.71) 



/ |/(x)|<ix < max \f(x)\(b — a). 

Jb xe[a,b] 



(3.72) 



Let I2 = J f(Qd£, noting that /(£) > 0. For any e > there exists some p large 
enough that /(£) < e/f for all < £ < £ Thus 





< |c 
- r 

< e. 



(3.73) 



I2 therefore also tends to zero in the limit p — > 00. Thus in the limit p — > oo, the 
expectation value ^ ^^ y — >■ 0, proving Lemma 13.4.1 1 
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Proof of Lemma I3A21 

If we take the forgetting time tj — > and keep the discreteness scale d p i fixed then / 



oo. 



First note that for all £, cosh£ > sinh£. The expectation value, Equation 13.671 can 
therefore be written 
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(3.74) 



Thus in the limit r/ — >■ 0, the expectation value — >■ oo proving Lemma 13.4.21 
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3.5 Consequences and bounds 

It is important to reiterate that the swerves diffusion equation is very general: any 
Lorentz invariant process for massive particles that results in fluctuations in energy- 
momentum will be governed by this equation. Deriving such an equation is of lit- 
tle use, however, unless we investigate the observable consequences. Ideally, a new 
phenomenological model should agree with all current observational data and, in ad- 
dition, make predictions for new, feasible, experiments. In reality, models will al- 
ways contain free parameters that limit their predictive powers. The swerves diffu- 
sion equation contains only one free parameter, the diffusion constant k, and is thus 
a strong phenomenological model. The usefulness of this model is unfortunately still 
limited. It will never be possible to prove that swerves do not occur - null obser- 
vations only allow the parameter k to be constrained to smaller values. It is hoped, 
of course, that experiments will actually observe the swerves effect; or that swerves 
will provide an explanation for some already observed, unexplained phenomena. The 
difficulty is (as discussed in the Introduction) that quantum gravitational effects are 
inherently very small and to make positive observations we must think of ways in 
which these tiny effects could have been amplified. For swerves the first step is thus 
to consider how the magnitude of the diffusion constant can be constrained by current 
data. Strong constraints have been placed on k by both Dowker et al. |3] and Kaloper 
and Mattingly ||30| . These constraints are reviewed in this section. Lacking the in- 
spiration of an astrophysical or cosmological problem swerves may solve, it seems 
unnecessary to seek tighter bounds at this stage. 

For simplicity and consistency with 130| I will work in terms of a scalar dis- 
tribution p = pt/y/g rather than the scalar density pt- In terms of the scalar p the 
homogeneous cosmic time diffusion equation becomes: 



The swerves model predicts the spontaneous heating of a nonrelativistic gas. To 
see this, consider the nonrelativistic limit of Equation 13.751 Working in polar coordi- 




(3.75) 
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nates on H 3 



and thus 
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In the nonrelativistic limit p « m this becomes 



dp 
dt 



p 2 sin I 



d 

— I p~sm„ 
ap Y 



dfl V d9 



]_0_ ( 2 dp 
p 2 dp Y <9p 

i a 2 p 



+ 



dp 

) 

d ( 1 dp 



+ 



smtf c 

1 d ( s[n6 dp 
p 2 sin 6 d6 \ d6 



+ 



p 2 sin 2 i 



kV 2 p, 



(3.77) 



(3.78) 



where V 2 is the standard Laplacian on M 3 , i.e. in the nonrelativistic limit the swerves 
diffusion equation reduces to the standard three dimensional diffusion equation. The 
general causal solution to Equation 13 . 78 1 can be derived using Greens's function tech- 
niques. The result is 



(Airkt) 
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4kt 



/(pVp' , 



(3.79) 



where /(p) is the initial distribution /(p) = p(p, 0). Consider an initial distribution of 
particles at rest in our frame, i.e. /(p) = <5 3 (p). In this case, the solution of the diffusion 
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equation simplifies to 



1 (P"P n2 



^ = I ^^ eXP l""J' 3(P ' )dV 



= (4vrA:t)- 3 / 2 exp^-^j . (3.80) 

This is just the Maxwell-Boltzmann distribution with a temperature T = where 
is the Boltzmann constant. The temperature is dependent on the length of time 
the diffusion process has been occurring - in other words, the nonrelativistic limit of 
the swerves diffusion equation describes a spontaneous heating of a thermal nonrel- 
ativistic gas. It is also clear that an initially thermal distribution will remain thermal. 
Suppose the initial distribution is a thermal distribution with a temperature To: 



/(p) = (2vrmfc i? To)- J ^exp - r - . (3.81) 
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where T = To + 2kt/mkB, i.e. a thermal distribution remains thermal, with the tem- 
perature changing at a rate of dT/dt = 2k /inks- 

Dowker et al. |3] placed a bound on the diffusion constant k by considering the 
spontaneous heating (or rather, the lack of spontaneous heating) of hydrogen gas in 
the laboratory. They assumed that a heating rate of 10~ 6 degrees a second would have 
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already been detected if it existed. This resulted in the bound 

rnks dT 



k 



2 dt 

< KT 56 % 2 mV 3 , (3.83) 



or equivalently 

k < lO~ u GeV 3 . (3.84) 
The average energy of a thermal distribution at temperature T is given by 

(E) = (3.85) 

Constraints on k can thus also be expressed in terms of the maximum energy gain 
allowed 

mA(E) 

k < (3.86) 

3At v ' 

To more tightly constrain k it is therefore necessary to consider systems that have 
been around for a long time without much gain in energy. Kaloper and Mattingly [30J 
looked at astrophysical clouds that have existed for roughly the age of the Milky Way 
without significant energy gain. They found, however, that energy loss due to radi- 
ation cannot be neglected and the result is a bound on k roughly comparable to that 
given above. They obtained a much stronger bound of k < 10 _61 GeV 3 by considering 
cosmic neutrinos. When heated by swerves, the cosmic neutrino background would 
appear as hot dark matter. The diffusion constant is thus bounded by observations 
that suggest dark matter is mainly cold. 

In attempting to bound the diffusion parameter, k, it is assumed that k is the same 
for all particles. Other than for simplicity, there is no reason to believe this is the case. 
A strong bound on k for neutrinos does not, therefore, rule out observable swerves for 
other particles. Without an underlying model for particles on a causal set (a model for 
the particles themselves, not just their propagation) it is unlikely any progress will be 
made determining the dependencies of k. 

It is also interesting to note that, although Lorentz invariant, swerves could give 
rise to an apparent signal of Lorentz violation 1135 1 . This is because many investiga- 
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tions into Lorentz violation rely not only on the details of the Lorentz violating model, 
but also on the assumption of energy-momentum conservation. A violation of energy- 
momentum conservation could be incorrectly interpreted as a signal of Lorentz viola- 
tion. Consider, for example, a time-of-flight experiment where a particle is measured 
to arrive with energy E. If the particle is subject to swerves, the arrival energy, E, is 
not the energy of the particle throughout its propagation, and the propagation time 
will therefore be miscalculated if swerves are neglected. Mattingly |35| concludes, 
however, that any false signal of Lorentz violation that appears due to the energy- 
momentum violation of swerves would be considerably smaller than current limits of 
experimental sensitivity. 
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3.6 Simulations of particle models 

The diffusion equations derived in this chapter do not rely on any particular model 
for particle propagation on a causal set. The diffusion parameter k is, however, likely 
to depend on the underlying propagation mechanism. Simulating the particle models 
discussed in Section 13.11 shows explicitly that the models lead to diffusion like be- 
haviour. It also allows an initial investigation into how the continuum diffusion pa- 
rameter might depend on characteristics of the underlying physical process. 

To simulate the particle models causal sets are constructed via sprinkling into 
Minkowski spacetime (as discussed in Section 12.31) , using the Cactus numerical rela- 
tivity code |36| and the CausalSets arrangement written by David Rideout. In this en- 
vironment, code was developed to implement the three particle trajectory algorithms 
discussed in Section 13.11 The causal sets generated in this manner can of course be 
considered as partially ordered sets independent of Minkowski spacetime. Here the 
information about the coordinates of the points in Minkowski spacetime will be re- 
tained to allow the causal sets to be easily visualised. The particular frame to which 
the coordinates of the sprinkled points in Minkowski spacetime refer will be called the 
embedding frame. 

The results discussed here are, for the most part, in 1+1 dimensions. The code is, 
however, completely general and can also be run in 3+1 dimensions (example trajecto- 
ries in 3+1 for each particle model are shown in Figures l331 and l3lj)) . While simulations 
in 3+1 are more physically relevant, they require much larger causal sets to obtain the 
same discreteness scale and the necessary time and computer resources were unfortu- 
nately not available during the writing of this thesis. Simulations in 1+1 dimensions 
in fact prove to be more than adequate to investigate the link between the continuum 
and discrete processes. 

Several important aspects of the simulations need to be noted. Firstly, an element 
is added to the sprinkling at the origin to allow a fixed beginning point for the trajec- 
tories. Also, the models all assume a partial particle trajectory already exists and is 
being extended - in these simulations the trajectories start with a single element and 
thus the development of the first part of the trajectory needs to be carefully considered. 
For Model 1, the swerves model, the matter is very simple. An initial momentum, p$ 
is defined and no further information is needed to construct the future trajectory. In 
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Figure 3.5: An example trajectory using Model 1 in 3+1 dimensions (one spatial dimension 
is suppressed). The causal set here has size N = 32768. The trajectory is constructed using 
Tf = 0.37 in embedding coordinates. 

these simulations po is taken as po = 0, i.e. the particle is initially at rest in the embed- 
ding frame. 

Model 2 relies on knowledge of e ra _i as well as e n to determine e n+ i. Unfortu- 
nately it is not possible to just define an e_i to begin the trajectory. In this case, when 
constructing the first step, the condition d(e n _i, e n +i) < 2n/ is simply ignored. The 
element e\ is chosen such that it maximises d(eo, e\) < nt. Points satisfying this con- 
dition will be close to the hyperbola defined by a proper time r = \ft 2 — x 2 = ntd p i, 
where d p i is the discreteness scale of the causal set. Again, the particle should be ini- 
tially at rest (or as close as possible) in the embedding frame. To impose this, instead 
of choosing e\ at random from those eligible, the element with the smallest time coor- 
dinate in the embedding is chosen. The rest of the trajectory is constructed as per the 
algorithm. 

Model 3 presents a slightly more complicated situation. Here the future trajectory 
from eo depends on elements e_ n/ , . . . , e_i. In simulations of Model 3, 1 first construct 
a longest chain between the element at t = 0, x = and an element att = 0.25, x = 
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(a) Model 2, nj = 5. 




(b) Model 3, n f = 5. 



Figure 3.6: Example intrinsic trajectories in 3+1 dimensions (one spatial dimension is sup- 
pressed). The causal set here has size N = 32768. 
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(also artificially added to the sprinkling). This longest chain forms the beginning of the 
trajectory: essentially it imposes the condition that the beginning of the trajectory is as 
close to geodesic as possible, and again the particle is initially at rest in the embedding 
frame. The rest of the trajectory is constructed as per the algorithm. Note that the 
initial longest chain may have less than n j elements, in which case minimisation is 
over all elements until nj elements are reached. 

A second important point to note about the simulations is scale: the causal set 
models give rise to the continuum diffusion equation when 1 < < n f < < I where / 
is the trajectory length. In other words, the causal set must be large enough that the 
trajectory involves many nj steps, and nj itself must be much larger than 1. This is 
quite difficult to obtain. Suppose n f = 10 and we want the trajectory length to be 10n y 
steps (or in the case of Model 3, 100 elements). In 1+1 dimensions this requires a causal 
set of roughly 10000 elements - quite manageable. If nj = 100, however, a causal set 
of roughly 1 x 10 8 is needed for I = lOOnjI^ This is certainly not manageable with the 
current code, so it is not really possible to explore realistic values of n f and trajectory 
length. Despite this limitation the simulations clearly show diffusion behaviour, as 
will be demonstrated below. 

3.6.1 Example trajectories in 1+1 dimensions 

Simply constructing an example trajectory for each of the three models in 1+1 dimen- 
sions reveals something of the relationship between the parameters of the model and 
the diffusion parameter in the continuum. Figure 13.71 shows a trajectory constructed 
in a causal set of size N = 32768, using the swerves momentum model, Model l|^| 
For this trajectory -ry = 0.1105 in embedding coordinates, and the final length of the 
trajectory is Z = 16 elements. Clearly this trajectory does not have a forgetting time, 
Tf, many orders of magnitude greater than the discreteness scale, or a total trajectory 
length much greater than Tf. The result is, however, very very close to the expected 
straight line in the continuum. 

Model 2 also gives a reasonably straight line for comparable parameters. Figure l3T8l 
shows a trajectory constructed in a causal set with N = 32768 and nj = 20. Note that 

2 In 3+1 dimensions you would need the absurd number of 1 x 10 16 elements. 

3 The causal sets are generated by a Poisson process with a mean of N, the actual number of elements 
varies. 
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Figure 3.7: Example trajectory using Model 1 in 1+1 dimensions: N = 32768, 17 = 0.1105, 
I = 16. 
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Figure 3.8: Example trajectory using Model 2 in 1+1 dimensions: N = 32768, n/ = 20, 1 = 13. 
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-0.5 -0.4 -0.3 -0.2 -0.1 0.1 0.2 0.3 0.4 0.5 
(b) Model 3, n f = 40,1 = 154. 

Figure 3.9: Example Model 3 trajectories in 1+1 dimensions: N = 32768. 
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Figure 3.10: 30 Model 3 trajectories in 1+1 dimensions: N = 32768, n } = 40. 

nj = 20 is approximately equivalent to tj = 0.1105 for a causal set of size N = 32768. 
Recall, however, that Model 2 actually relies on information 2n j into the past of the 
trajectory: it would be more correct to compare a nj = 10 Model 2 trajectory with a 
tj = 0.1105 Model 1 trajectory. Regardless, it is clear that with roughly comparable 
forgetting times Model 2 results in more fluctuations in position and momentum. 

A comparison of the two intrinsic models, Model 2 and Model 3, is also interesting. 
In Figure [3T91 two trajectories constructed using Model 3 are shown: one with nf = 20 
and one with nj = 40. As mentioned above, Model 2 relies on information 2ny steps 
into the past, so we might expect similar trajectories from n f = 20 Model 2 and n f = 40 
Model 3. The trajectories from Model 3, however, swerve much more. Given that the 
values of n f and I are not really in the realm 1 << n j << /, where the continuum 
diffusion is expected to hold, care should be taken in drawing conclusions about the 
continuum behaviour of Models 2 and 3. Even so, it seems reasonable to conclude 
that a simple universal relation between the forgetting number of an intrinsic particle 
model and the diffusion parameter in the continuum does not exist (a concrete relation 
between the forgetting time of Model 1 and the diffusion parameter will be found in 
Section |3.6.2|) . Of course, the models here are not physically realistic. It seems likely, 
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however, that when physically realistic models are developed, the determination of 
the diffusion constant k from the 'fundamental' parameters of the underlying model 
on the causal set, will be model dependent. 

The diffusion-like behaviour of the models is most obvious if multiple trajectories 
are generated. Recall that the intrinsic models do not define unique trajectories - an 
element is chosen at random from the eligible elements at each step. For a fixed causal 
set and initial trajectory, multiple trajectories can thus be generated. Figure [37TOl shows 
30 (distinct) Model 3 trajectories in a fixed causal set, N = 32768, nj = 40. Note that 
from t = to t = 0.25 the trajectory is a fixed longest chain for all 30 trajectories. 
Although diffusion-like behaviour is clear, it is difficult to determine whether it is 
exactly the diffusion given by the swerves diffusion equation. The values of N and n j 
are nowhere close to the continuum limit, and 30 trajectories is simply insufficient. To 
compare more directly with the swerves diffusion equation, I will look instead at the 
non-intrinsic model, Model 1. 

3.6.2 A close investigation of Model 1 
Simulating trajectories 

Unlike the intrinsic models, for a given causal set and fixed initial position and mo- 
mentum Model 1 defines a unique trajectory. To say that Model 1 results in diffusion 
is, in a sense, saying that the exact underlying causal set is unknown. To simulate this 
I generate many different sprinklings into the same region of Minkowski spacetime, 
and calculate the unique trajectory in each. As above, the starting point of each tra- 
jectory is fixed at x = 0, t = 0, an extra element artificially added to each sprinkling. 
For each trajectory a final position and momentum can be calculated at tf ina i = 0.95. 
This is done by identifying elements e n and e n+ \ in the trajectory such that e n is the 
last element with t < t fi na i and e n+ i is the first element with t > t fi na i . Note that there 
is a vanishingly small probability that any element has t = tn na i. The final position 
is found by linear interpolation between e n and e n+ \. The final momentum is simply 
the usual momentum at e n+ i, i.e. it is proportional to the vector between e n and e n+ \. 

The first step in the analysis is to reject some trajectories - if a trajectory is close 
to the boundary of the region of Minkowski spacetime at any point, it will 'bounce' 
back and distort the results. An example of such a trajectory is given in Figure 13.111 
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Figure 3.11: An invalid Model 1 trajectory, N = 32768, 77 = 0.02. 

Provided only a small fraction of the trajectories are invalid, the invalid ones can be 
safely rejected without much loss of information about the trajectory distribution. The 
parameter values investigated in this section are chosen such that there are indeed 
few, if any, invalid trajectories. 

Once the invalid trajectories have been removed, histograms of final position and 
final momentum can be generated from the remaining trajectories. These histograms 
clearly show diffusion behaviour, but to truly test the model the results must be com- 
pared with the distributions expected from the swerves diffusion equation. 

The swerves equation 

The inhomogeneous, 1+1-dimensional swerves diffusion equation is 

dpt _ ~P dpt | h & i ^m 2 + p 2 dp t \ ^ 3g7 ^ 
dt sjm 2 + p 2 dx dp I m dp J 

As this does not have an exact solution, it must be evolved numerically. For the results 
shown here this numerical evolution is conducted with NDSolve in Mathematica. 
The initial distribution of the simulated trajectories is a ^-function: all trajectories be- 
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gin with x = and p = 0. To numerically evolve Equation 13.871 with a 5-function 
initial distribution is not possible - instead, the initial distribution is taken to be a 
highly peaked Gaussian 



imation to a histogram of the initial position-momentum data constructed with the 
same bin size as used for the final data, and such that Equation l3.87l evolves smoothly 
without numerical artifacts. To obtain final position and momentum distributions the 
evolved pt(%,P,t final) is summed over momentum and position, respectively. 

The matter of units 

To compare the simulated trajectories with the swerves diffusion equation a consis- 
tent set of units needs to be chosen: Equation 13.871 above, is expressed in terms of 
Planck units; the simulated trajectories are expressed in terms of 'embedding' coordi- 
nates. As mentioned earlier in this thesis, the discreteness scale in causal set theory 
is expected to be of the order of the Planck scale. For concreteness I assume the two 
are the same, and thus Equation 13.871 can be considered as expressed in 'discreteness 
units'. A discreteness length d p i = y/V/N, where V is the volume of the region of 
Minkowski spacetime in embedding units and N is the mean causal set size, can be 
calculated from the simulated causal sets. The position, time, momentum, and mass 
in the simulations are then re-expressed in terms of discreteness units. 

Comparing the simulations to the swerves equation 

Equation 13.871 is evolved for a range of values of k and the resulting distributions are 
compared to the position and momentum histograms for given values of the model 
parameters (r/, N, and particle mass to). A best fit value of k for each set of model 
parameters is found by minimising the reduced x 2 '- 




(3.88) 



Here, o 2 x and a 2 are not fixed: they are chosen such that Equation l3.88l is a good approx- 




(3.89) 



Xred X 7/: 



(3.90) 
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where Oj is the observed frequency, Ei is the expected frequency (i.e. that given by 
the evolution of the diffusion equation), and / is the number of degrees of freedom 
(here, / = number of data points - 1). The reduced x 2 is calculated for the distribution 
in momentum rather than the full momentum-position distribution. The momentum 
distribution is chosen for comparison as it drives the position diffusion. Care must 
be taken when calculating the value of the reduced x 2 ~ ^ a significant proportion 
of the frequency values are less than five, x 2 is no longer a good measure of fit (see, 
for example 1371). In situations where the frequency in a given bin is less than five, 
multiple bins are combined to resolve this problem. A reduced x 2 value of x 2 e d ~ 
is usually taken to indicate a good fit (see, for example 138]), i.e. the evolved diffusion 
equation is a good description of the simulated trajectories. 

Results 

The goal of these simulations is not just to demonstrate that the model gives diffusion 
behaviour, but also to determine the relationship between the continuum diffusion 
parameter k and the parameters of the underlying discrete model. The relationship 
between k and the forgetting time Tt can be determined by running simulations for a 
range of values of 77 (with fixed iV and m) and determining the best fit value of k in 
each case. 

A causal set of size N = 32768 was used in these simulations. The particle mass 
in embedding units is taken to be m = 1. Five hundred trajectories were evolved for 
each value of 77 in 77 = {0.03,0.035,0.04,0.045,0.05,0.055,0.06,0.07,0.08,0.09,0.1} 
(in embedding units). The results are summarised in Table I37T1 Several things need to 
be considered when interpreting these results, as discussed below. 

Table I3T1 shows the simulated values of 77 in both embedding units and multiples 
of the discreteness length d v %. In units of d v \, 77 varies between 5.4 and 18 - clearly 
not the many orders of magnitude greater than d p i that is expected in reality. The 
average trajectory length (this is the number of 'steps' in the trajectory rather than 
the number of elements) is also well out of the realistic range, reaching a maximum 
of only 48. Unfortunately realistic values are simply not feasible with the current 
computing facilities. What is remarkable is that, even with these values of 77 and 
N, the distribution of the simulated trajectories is extremely well modelled by the 
diffusion equation. The best fit values of the diffusion parameter and corresponding 
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y (d p i) 


Number of invalid 
trrHprtorips 


Averatre traiectorv 
length 


best fit k 


~X-red 


0.03 


5.4 


10 


48 


1.2 x 10~ 8 


1.5 


0.035 


6.3 


2 


42 


5.5 x 10~ 9 


1.0 


0.04 


7.2 





37 


2.8 x 10~ 9 


0.59 


0.045 


8.1 





33 


1.5 x 10~ 9 


0.64 


0.05 


9.1 





30 


9.2 x 10~ 10 


1.0 


0.055 


10 





27 


5.0 x 10~ 10 


2.4 


0.06 


11 





25 


3.7 x 10~ 10 


1.6 


0.07 


13 





22 


1.7 x 10~ 10 


0.56 


0.08 


14 





19 


9.9 x 10~ n 


0.55 


0.09 


16 





17 


4.7 x 10~ n 


2.0 


0.1 


18 





16 


2.6 x 10~ n 


2.4 



Table 3.1: Analysis of Model 1 trajectories for various values of Tf, with fixed N and m. 



xfed values are given in Table 13.1 1 As mentioned earlier, x 2 re d ~ 1 is usually considered 
a good fit. It is easy to see just how good the fit is when the evolved diffusion equation 
is plotted together with the simulation histograms. Examples of both the position and 
momentum histograms for four of the values of Tf are shown in Figure 13.121 The 
particular values of Tf shown are not special in any way, and are merely chosen to 
illustrate the results. 

As Table l3Tl shows, invalid trajectories are an issue only for the smallest two values 
of Tf and even in the worst case only comprise 2% of the 500 trajectories. As such, they 
don't influence the results given here. They do, however, become significant for values 
of Tf smaller than 0.03 lessening the utility of simulations for those values. 

Although the best fit values of k and corresponding Xred calculated here are fairly 
robust, it should be noted that the best fit values of k are determined by calculating 
xf e d f° r various values of k and choosing the one that minimises x 2 re d- A truly system- 
atic search of the parameter space is not undertaken, and a measure of the error in k is 
not calculated. The exact values of xfed d e P en d slightly on the bin width chosen when 
generating the histograms. They also vary slightly if a x and a p in the initial Gaussian 
distribution are changed. A further point is that for some Tf there are a range of values 
of k that give x\ e d ~ 1 an d thus all could be considered a good fit. For the results that 
I have obtained, k vs. Tf is plotted in Figure [37l3l The dependence of k on rj is clear in 
the In — In plot, where the best fit line to the data, In k = —5.02 In Tf — 9.8 is also shown. 
Despite the considerations mentioned above, it seems clear that k ~ 1/t5 . 
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(a) Tf = 0.04. 
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(b) T f = 0.05. 
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(c) t s = 0.06. 
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(d) T S =0.08. 

Figure 3.12: Example histograms with best fit numerical solutions. The position distribution 
is shown on the left, momentum on the right. 
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(b) In k vs. In 77 with the best fit line -5.02 In r/ - 9.8. 
Figure 3.13: Comparison of the diffusion parameter k and forgetting time Tf, for fixed TV and 



m. 
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N 


dpi 


Number of invalid 
trajectories 


Average trajectory 
length 


best fit k 




4096 


0.016 


12 


19 


0.034 


0.64 


8192 


0.011 





19 


0.0088 


1.0 


16384 


0.0078 





19 


0.0021 


0.45 


32768 


0.0055 





19 


0.00056 


0.55 



Table 3.2: Analysis of Model 1 trajectories for various values of d p i (in embedding units). 

The diffusion parameter k also depends on the particle mass m. In Model 1 the 
trajectory is independent of the particle mass - the mass is only used to correctly 
normalise the momentum. The momentum histogram for trajectories of a particle 
of mass m' = am can therefore be obtained simply by scaling the original momentum 
histogram by a. Looking at Equation 13.871 if both m and p are scaled by a and k is 
scaled by a 2 the equation remains unchanged. Thus it can be determined, without 
reanalysing the simulations, that k ~ m 2 . 

There is one remaining parameter that k can depend on: the discreteness length 
itself. From dimensional analysis of Equation l3.87l it can be deduced that k has units of 
L~ 3 (since Planck units are used, L = T = M" 1 ). Since k ~ m 2 Tj 5 the dependence on 
dpi must be k ~ dt v This dependence can be checked through simulations. In this case 
everything needs to be expressed in terms of embedding units rather than discreteness 
units. Five hundred trajectories are calculated for four different values of N (and thus 
dpi) and the same values of forgetting time Tf = 0.08 and mass m = 1. The value 
of Tf was chosen to avoid too many invalid trajectories in the causal sets considered. 
As above, the best fit values of k and the corresponding Xred can " 3e calculated. The 
results are summarised in Table 13.21 Again, a In — In plot (see Figure 3.14(a)| ) allows 



the dependence of k on d p i to be determined. The result is as expected: k ~ d. 



pi- 



The analysis of the simulations leads to the relation k ~ m 2 d^ l /Tj. The constant of 



proportionality can in fact also be determined. Looking back at the results for varying 
Tf and noting that in terms of discreteness units m = 0.0055 and d p i = 1 the relation 
k vs. m 2 /r| can be plotted (see Figure 3.14(b)| ). The best fit line is found to have the 



equation k = 1.9 m /rf — 2.5 x 10 . Thus, it can be concluded that 



k « 2m 2 4/r / 5 



(3.91) 



Causal set theory contains not only a discreteness scale, but also a hypothesised 
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Figure 3.14: The relationship between the diffusion parameter k and underlying model pa- 
rameters d p i, to, and r/. 
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'nonlocality' scale below which physics cannot remain local |39l l40l . The forgetting 
parameter that appears in the massive particle models may be a measure of this non- 
locality scale. Sorkin |39| estimates the nonlocality scale to be of the order 10 20 in 
Planck units. Although the simulations here use forgetting times nowhere near this 
estimate, and there is no theoretical reason why the forgetting time should be exactly 
the nonlocality scale, it is interesting to consider the value of the diffusion constant 
that arises if the above results are extrapolated. Taking the particle to be a proton with 
mass m = 1CT 20 , r f = 10 20 , and d pl = 1, Equation l3T9Tl leads to the estimate k « 1CT 140 . 

As discussed in Section l331 tight constraints have already been placed on k through 
observations. In particular (see Equation l3.83|) k < 10 _56 %m 2 s -3 ~ 1CT 102 due to the 
lack of spontaneous heating of hydrogen in the laboratory. The diffusion parameter 
estimated from the nonlocality scale is thus many orders of magnitude smaller than 
the current constraint (the factor of two mass difference between the proton and a hy- 
drogen molecule is inconsequential here). This rough calculation should not be given 
too much weight, but it does suggest that k could be well outside our current observ- 
able range. 
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3.7 Concluding remarks 

Following the proposal by Dowker et al. of the swerves model, and their derivation of 
the proper time swerves diffusion equation, the question arose of how an underlying 
model of particle behaviour on a causal set could be formulated without relying on 
information about the approximating continuum spacetime. The association between 
longest chains in a causal set and geodesies in the continuum allowed such models to 
be created, two examples were given in Section [37T1 The intrinsic models described in 
this thesis are by no means the only options, and simply give an idea of the wealth 
of possibilities available. The models proposed here (including the original swerves 
model) make no attempt to be physically realistic. A true description of particle mo- 
tion in causal sets will have to await many further developments in the theory. The 
first step might involve developing a quantum mechanical rule for the trajectory evo- 
lution, or modelling particles as wavepackets. The models proposed do capture the 
crucial features of causal set theory - they could, in a sense, be interpreted as a coarse 
grained description. Just as modelling a gas as an unrealistic collection of hard spheres 
proved enormously useful, so too do these models prove useful in determining con- 
tinuum behaviour. 

The proper time swerves diffusion equation, first given in [3| and derived in full in 
Section 13721 is the unique continuum diffusion equation given an underlying Marko- 
vian, Poincare invariant process causing random fluctuations in momentum. It is in- 
dependent of the particular underlying models discussed - those models and many 
more would result in continuum behaviour described by the swerves diffusion equa- 
tion. One query that remained, was whether a given model was guaranteed to result 
in a finite diffusion constant. In Section [374] it was shown analytically that the swerves 
model does lead to a finite diffusion constant. 

The swerves diffusion equation, expressed in terms of proper time, is difficult to 
compare with observations. Although bounds had been placed on the diffusion pa- 
rameter by considering the non-relativistic limit, bounds in the fully relativistic case 
could not formerly be considered. The full inhomogeneous swerves diffusion equa- 
tion in terms of cosmic time was derived in Section 13731 With this equation complete 
comparison with observations will be possible. 

After exploring swerves analytically, the second task of this chapter was to explore 
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the models numerically. The underlying models give rise to diffusion in the contin- 
uum limit. With simulations limited to quite small causal sets (far from the infinite 
sprinkling density of the continuum limit) it was not expected that the full diffusion 
behaviour of the models would be apparent. The results were, in fact, far better than 
anticipated. All three models clearly show diffusion like behaviour. A close investiga- 
tion of the swerves model shows the model behaviour perfectly matches the swerves 
diffusion equation. It was even possible to determine an expression for the continuum 
diffusion parameter in terms of the parameters of the underlying causal set model. 
That the diffusion behaviour occurs even this far from the continuum limit is very 
interesting. There are also wider implications for causal set simulations - the results 
here show that despite computational limits on causal set sizes, useful information 
about the continuum appearance of discrete spacetime can be obtained. 



Chapter 4 



Massless particles 



Die ganzen 50 Jahre bewusster Griibelei haben mich der Antwort der Frage 
„Was sind Lickhtquanten" nicht naher gebracht. Heute glaubt zwar jeder 
Lump, er wisse es, aber er tauscht sich. 

A total of fifty years of careful reflection has not brought me any closer to 
answering the question 'what are the quanta of light?'. It's true, today, that 
many a fool believes he knows the answer, but he's mistaken. 

Einstein to Besso, 12th December 19510 

An underlying spacetime discreteness results in diffusion in momentum and space- 
time for massive particles. Does a similar diffusion occur for massless particles? For 
massive particles concrete models of particle propagation on a causal set were con- 
structed, motivating the derivation of the swerves diffusion equation. For massless 
particles such underlying models are unfortunately elusive. Consider a sprinkling 
into Minkowski spacetime: for any given element, p, there will almost surely be no 
element sprinkled on the future lightcone of p. A trajectory constructed, as in the mas- 
sive case, by jumping from one element p to another q, where p -< q, will rarely be 
close to null. A close analogue to the future lightcone of p can, however, be defined 
simply as the set of all elements q linked to p: p -<*q. These elements q lie near the 
hyperboloid with proper time one Planck unit from p. This suggests that the simplest 
analogue of a null ray is a single link. More generally, p -< q is a good approximation to 
a section of a null ray if the causal interval between p and q is totally ordered |40l[T2| , 
but it is still hard to see how a discrete Markovian process resulting in a close-to-null 
trajectory might be defined. It may be that options such as allowing particles to also 

Published in |41 1, English translation from the French translation in FiTI courtesy of Dr. R. M. Pollard. 
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move backward in time need to be considered. More generally, it is not surprising that 
attempts to model massless particles as classical point particles run into trouble. It is 
hoped that the study of massless fields on a causal set | 39l |42| and general quantum 
field theory on a causal set |32ll33 ] will eventually lead to a model for massless particle 
propagation at a discrete level. 

The current lack of an underlying model does not, fortunately, exclude the inves- 
tigation of the potential effect of discreteness on photons in the continuum approxi- 
mation. A massless diffusion equation can be derived in one of two ways: following 
the stochastic evolution on a manifold of states procedure, as for the massive case; 
or, simply taking a m — > limit of the massive particle diffusion equation. As will 
be discussed in the following section, it turns out that the second method gives an 
incomplete result. 

In Section 14.11 the massless particle diffusion equation in terms of affine time is 
derived. As for massive particles, in order to make comparisons with observations an 
equation in terms of cosmic time is needed - such an equation in derived in Section l4~2l 

The resulting massless particle model has two free parameters. In Section l4~3l the 
effect of massless particle diffusion on the cosmic microwave background is studied, 
and bounds are placed on the values of the two free parameters. Finally, the conse- 
quences of the model for astrophysical spectra are discussed in Section 14.41 where it 
is found that the strong bounds from the cosmic microwave background rule out ob- 
servable effects. The work in Sections l4. 1 [ 14.21 l4~3l was conducted in collaboration with 
Fay Dowker and Rafael Sorkin, and appears in [29]. 

Note that the massless particles considered in this chapter are assumed to have no 
polarisation. The issue of polarisation is addressed in depth in Chapter [5] 
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4.1 Diffusion in affine time 

The state space in the massless case differs from the massive case. For massive parti- 
cles there is a probability distribution on M 4 x H 3 . For massless particles H 3 becomes 
the light cone in momentum space defined by = 0. This cone will be denoted Hq. 
Hereafter the momentum of a massless particle will be denoted Assuming that the 
photons under consideration are well described in a geometrical optics approxima- 
tion (and thus have definite spacetime worldlines and momenta), the state space will 
be M 4 x Hq. 

Since proper time vanishes along a lightlike worldline, it is clearly not a suitable 
time parameter for a massless diffusion process. Instead, an affine time, A, defined by 

dx^ = k^d\ , (4.1) 

is used to parametrise the photon worldline. 

In the massive particle case the density of microstates, n, was proportional to the 
determinant of the metric on the state space: n oc yj~g. In the massless case, the induced 
metric on Hq is degenerate. A volume element on Hq can, however, still be defined. 
The four dimensional volume element d 4 k of momentum space, together with the 
masslessness constraint, k^k^ = 0, gives the invariant volume element d 4 k 5{k^k tx ) = 
d s k/2k° on Hq, i.e. n oc l/k° in Cartesian coordinates. It turns out to be more useful, 
however, to work in polar coordinates on Hq: {k, 8, cf)} where k is the magnitude of the 
three momentum and 6 and <p are the usual polar angles in momentum space. In these 
coordinates, the density of states is n oc k sin 0. 

As with the massive particle case, the derivation of the massless diffusion equation 
begins with the current and continuity equations, Equation l3.5l and Equation l3.6l The 
first step is to determine K AB and u A . To begin, note that there is a unique, up to a 
constant factor, invariant vector field on Hq which is the momentum itself, k a , i.e. the 
vector with components (k, 0, 0) in polar coordinates. This is absent in the massive 
case, where the momentum vector does not lie in the mass shell. Also, although there 
is no invariant metric on Hq, there is an invariant symmetric two-tensor, k a k b , again 
unique up to a constant factor. 
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From Equation B.llt the spacetime component of K 



is 




lim k^k 



Ax^Ax 



2AA 



AA 



,v 



AA->0+ 



0. 



(4.2) 



K AB is positive semidefinite so K^ a = 0. Finally, the requirement that K AB be Lorentz 
invariant and translation invariant fixes the remaining component K ab to be propor- 
tional to k a k b : 



where d\ > is a constant. 

The vector u A is easily determined by looking individually at the components in 
spacetime and momentum space. As in the massive particle case, the spacetime com- 
ponent u 11 = v 11 by Equation 13.161 and = k^ by Equation 13.121 In contrast to the 
massive case, there can be nonzero components of u A in the momentum space direc- 
tions because the momentum itself is an invariant vector. The momentum direction 
components are thus given by u a = d,2k a , where d 2 is a constant. Working in polar 
coordinates the 'position' vector k a on the cone Hq is simply (k, 0, 0) where k 2 = kg 2 . 
Thusu A = {k°, k 1 , k 2 , k 3 , d 2 k, 0,0) onM 4 x Hjj. 

Substituting the forms for K AB and u A into Equation l3.15l gives the massless par- 
ticle affine time equation: 



Lorentz invariance constrains the diffusion process such that, although there is 
diffusion in the energy of the photon, the direction of the photon is unchanged and 
it continues to propagate at the speed of light. Note that there are two parameters, 
d\ and d 2 , in this phenomenological model, rather than the single parameter of the 
massive particle model. There is not only a diffusion term but an independent drift 




(4.3) 




(4.4) 
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term, arising from the existence of an invariant vector on Hq. It will be shown further 
below that this extra term leads to the existence of power law equilibrium solutions of 
the diffusion equation. 

It is not unreasonable to expect that a diffusion equation for massless particles 
could be obtained simply by taking the m — > limit of the massive particle diffusion 
equation, Equation l3.25l This, however, results in Equation l4.4l with d,2 = 0, as there is 
no invariant vector in the massive case. It is interesting to note that there is this discon- 
tinuity in behaviour between very highly relativistic massive particles, and massless 
particles. It suggests that massive and massless particles are fundamentally different 
objects and will need to be described in completely different ways in causal set theory. 
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4.2 Diffusion in cosmic time 

As with massive particles, the massless particle diffusion equation is more useful for 
comparison with observations when expressed in terms of cosmic time. The massless 
cosmic time diffusion equation can be obtained by following the same argument as in 
the massive case (Section l331) . It is assumed that the massless particle trajectory begins 
at A = and t = 0. The first step is to define a new larger state space incorporating the 
affine time: M 1 = M 4 x ijj x K. In this larger space, a new current component can be 
defined: 

j\t,x i ,k a ,X) = p x . (4.5) 

The t component of the current J can be expressed in terms of J A , using Equations 13.51 
and EES 

j\t,x\k a ,X) = kJ x , (4.6) 
^Px = (4.7) 

The remaining components of the current can now be expressed in terms of J*. 

J l (t,x\k a ,X) = ^-J\ (4.8) 
' k 

J a (t,x\k a 1 X) = -dik a k h ksiu9d b ( 7 / - ) 
v ' \k z sm J 

+^d 2 k a , (4.9) 
k 

where Equation l3.31l has been used. In polar coordinates the vector k a = (k, 0, 0), and 
thus there is only one nonzero component of J k<l in the radial (energy) direction: 

(9 P 

J k (t,x\k a ,X) = -dxk — + (Mi + d 2 ) J 1 . (4.10) 

' ok 

Since the affine time along the particle trajectories is unobservable, it can be integrated 
over, and a new integrated current, J, can be defined. Integrating over the t compo- 
nent gives 

/•oo 

J\t,x\k a ) = / J t (t,x i ,k a ,X)dX. (4.11) 
J o 
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The remaining current components can be integrated and re-expressed in terms of J 1 : 

J\t,x\k a ) = jj\ (4.12) 
J k (t,x\k a ) = -d l k— + {2d l + d 2 )J t . (4.13) 

In the extended space M! , the continuity equation becomes 

d a J a = 0, (4.14) 

where X a = {X A ,X} are coordinates on M x R. Integrated over A the continuity 
equation becomes 

[J A ]§° + dj 1 + d l J i + d a J a = . (4.15) 

J x | a=o is zero for all i > and J A tends to zero as A approaches infinity for finite t. So 
for all t > 

d t J t + d i J l + d a J a = Q. (4.16) 
If J 1 is renamed p t the massless particle cosmic time equation can be written: 

% = -diP-dal* 

at 

= -J d m-(dx + d 2 )^ + dl k^. (4.17) 
4.2.1 Power law equilibrium solutions 

It is interesting to note that negative values of the drift parameter d 2 allow for power 
law equilibrium solutions of the massless particle cosmic time diffusion equation. For 
homogeneous distributions Equation l4.17l becomes 

-^ = -^l + d 2 ) — + d,k^. (4.18) 

Equilibrium solutions of this equation must satisfy 

° = -(*+*)^+**If- (419) 



:4.2 Diffusion in cosmic time 



78 



Equation l4.19l has a power law solution 

2d 1 +d 2 

pt oc . (4.20) 

If a small energy cut off k m i n is assumed, this solution is normalisable when the ex- 
ponent is less than -1. Since the diffusion parameter d\ is positive this can only occur 
for negative values of d%. It can be conjectured that if (2d\ + c?2)/di < —1 (equiva- 
lently je^l > 3di) any normalised distribution will tend to this power law equilibrium 
solution at late times. 
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4.3 Diffusion and drift in the cosmic microwave background 

In developing a phenomenological model the aim is to provide a model for currently 
unexplained observations or to suggest new observations that might be made to test 
a theory. Before proposing new observations or exotic experiments, however, one 
should constrain the model as tightly as possible based on what is already known. 
The massless particle momentum diffusion discussed in the preceding sections has 
two free parameters: a positive diffusion constant d\ and a 'drift' constant di, which 
may be either positive or negative. If the photons we encounter in our everyday lives 
experienced large amounts of diffusion or drift, it would of course have already been 
noticed. The motivation for the massless momentum diffusion model, however, lies 
in the underlying discreteness of causal set quantum gravity - the resulting effects are 
expected to only be apparent on very small scales, outside the range of current exper- 
iments. How then can the parameters be constrained? The key is to look at photons 
that have been travelling long enough for the diffusion and drift to 'accumulate' to an 
observable level. 

The cosmic microwave background (CMB) seems an ideal testing ground: not only 
are its photons the 'oldest' we can observe, but its spectrum has been determined with 
great precision. The photons in the CMB have been 'free streaming' since the surface 
of last scattering at a redshift of z ~ 1100 - approximately 13.7 billion years. When 
the universe became transparent at recombination, the photons that are now observed 
as the CMB would have had a blackbody spectrum with a temperature of 3000 K (see 
for example |43J). As the universe expanded the photons were both stretched and 
diluted, a process that neatly preserves the blackbody nature of the spectrum and 
simply lowers the temperature. Current observations of the CMB yield a temperature 
of 2.728 db 0.004 K and measure the spectrum to be Planckian (blackbody) over the 
2 — 21cm~ 1 frequency range to within a weighted root mean square (rms) deviation of 
only 50 parts per million (ppm) of the peak brightness l44l . Any diffusion and drift in 
energy would have distorted the CMB spectrum - the fact that the CMB photons have 
travelled so far but remained so perfectly thermal allows the parameters in the model 
to be constrained very tightly. 
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4.3.1 Modelling diffusion in the CMB 

The derivation of the massless cosmic time diffusion equation assumed spacetime to 
be Minkowskian. Although the expansion of the universe can be incorporated into 
the model (and will be in Section l4.3.2l below) a first bound on the parameters will be 
obtained by simply ignoring the expansion. In a non-expanding universe the temper- 
ature of the CMB would have remained constant from the surface of last scattering to 
today. To allow comparison with data, it will in fact be assumed that the CMB was 
emitted with the temperature observed today T = 2.72&K. 

If a blackbody spectrum is evolved according to the massless cosmic time diffu- 
sion equation, the final distribution will no longer be blackbody. The deviation from 
blackbody can be determined and, by requiring it to be within the current deviation 
allowed by observations, bounds can be placed on the drift and diffusion parameters. 
Specifically, the initial Planckian spectrum, expressed as a number density of photons 
per unit spatial volume per unit energy, is 

k 2 

p(k, t = 0) = 8^ -j— — , (4.21) 

exp (I) - 1 

with a temperature T = 2.72&K. This distribution is evolved over a time interval equal 
to that since the surface of last scattering (1 x 10 60 in Planck units) via the homogeneous 
massless cosmic time diffusion equation 

d Pt (A I A \ i J U^Pt (A m 

^ = - {dl + d2) M +dlk W^ (422) 
using the MAT LAB numerical pde solver pdepe. 

Boundary conditions 

Some questions must be addressed when choosing a boundary condition at k = for 
this integration. What happens to a photon as its momentum approaches zero? Do 
photons leak away through the tip of the null cone in momentum space? Physically, 
the photon concept employed by the model breaks down as the wavelength tends 
to infinity, because the geometrical optics approximation fails. This suggests that an 
'absorbing boundary condition', p(E) = 0, is appropriate at E = 0. As it happens, the 
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current is 



J = (2^ + d 2 )p t - dik—1- , 



(4.23) 



and thus (so long as dp t /dk remains finite) the absorbing boundary condition p t = is 
equivalent to the 'reflecting boundary condition', J = 0, at E = 0. In the simulations 
that follow, a reflecting boundary condition is used at the upper boundary k = k max . 

Calculating the deviation 

The evolved spectrum was converted from a number density per unit volume per unit 
frequency to a spectral radiance (energy per unit area per unit time per unit frequency 
per steradian) as used in the analysis of the COBE FIRAS data. This allowed the devi- 
ation from Planckian to be compared with the quoted 50ppm of the peak brightness. 
The first step was to fit a Planck spectrum to the evolved spectral radiance using the 
least squares method. Specifically, the MATLAB function fminsearch was used to 
find the temperature that minimised the sum of the squared differences between the 
evolved distribution and the fitted Planck spectrum. Looking for the best fit Planck 
spectrum rather than comparing with the initial 2.728K spectrum, allowed for the 
possibility that the diffusion changes the temperature of the CMB in a way that may 
be reconciled with observation. As it happens, it was found that the temperature of 
the best fit Planck spectrum was very close to the initial temperature in cases where 
the deviation is within the allowed tolerance. For example the choice of parameters 
di = 5 x 1CT 97 and d 2 = 1 x 1CT 96 gives a best fit temperature of 2.7281K, indistin- 
guishable from the current observed temperature of 2.728 ± 0. 004.fr. 

Finally, the rms deviation between the fitted Planckian spectrum and the evolved 
spectrum 



where s« is the value of the evolved spectrum at a point i, fa is the value of the fitted 
spectrum, and iV is the total number of data points, was calculated over the frequency 
range 2-21cmr 1 (energy range 4 x 10~ 23 -4 x 10 -22 J) with all points weighted equally. 
After calculating the peak brightness of the evolved spectrum, this rms deviation can 
be compared to the allowed tolerance of 50 parts per million of the peak brightness. 
This process was repeated for a range of values of the parameters d\ and d 2 . 




(4.24) 
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Results 

For the evolved spectrum to be close enough to Planckian to be within the allowed tol- 
erance, the deviation must be essentially invisible to the eye. Before placing bounds 
on d\ and cfo it is useful therefore to get a general feel for how a blackbody spectrum 
evolves under Equation l4.22l by looking at values of d\ and d,2 that result in deviations 
well outside the allowed tolerance. Figure |4~T1 shows the evolved distribution and best 
fit Planck spectrum for d\ = 1 x 1CT 93 , d2 = 1 x 1CT 93 . The positive drift parameter 
d% clearly shifts the spectrum to the right. With these large values of d\ and cfe the 
difference between the best fit Planck spectrum (2.851 K) and the initial Planck spec- 
trum (2.728 K) is quite considerable. However, no amount of temperature change will 
disguise the translation of the Planck spectrum - the diffusion and drift are clearly 
visible. 

To place bounds on the diffusion and drift constants they were first examined sep- 
arately, varying d\ with g?2 = and varying c?2 with d\ = 0. When d\ = Equation l4.22l 
can be solved exactly: 

pt(k,t) = po(k-(ht), (4.25) 

i.e. the spectrum just translates at a constant speed. For c?2 negative, this is inconsistent 
with the boundary condition p t = at k = 0. However, in this case one can implement 
an absorbing boundary condition trivially: simply cut off the translated distribution 
at k = 0. The exact solution is also useful as a method of checking the level of nu- 
merical error in the simulations. Since the deviations from Planckian that are being 
investigated are so small there is the possibility that the rms deviation calculated is 
significantly affected by numerical error. Fortunately, this is not the case - the differ- 
ence between the exact solution and the numerically evolved solution for d\ = is 
several orders of magnitude smaller than the deviation from Planckian^ For example, 
if efo = ±4 x 10~ 96 the rms deviation from the best fit Planck spectrum is 5 x 10~ 101 
(5 x 1CT 5 peak brightness) for both the exact and the numerical solution. The rms de- 
viation between the exact and numerical solutions is 4 x 10~ 104 . This also demonstrates 
that the p = boundary condition imposed on the numerical solution, although in- 



rerunning all the simulations in this section using the Mathematica numerical solver NDSolve 
yields identical results further suggesting that numerical error and integration step size do not have a 
significant effect on the bounds. 
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Figure 4.1: Comparison of evolved and best fit Planck spectrum for d\ = 1 x 10 93 , d 2 = 
1 x 10 ~ 93 . Energy k and spectral radiance are in Planck units. 
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Figure 4.2: The rms deviation of the simulated spectrum from Planckian as a proportion of the 
spectrum peak, varying d 2/ d\ = 0. 

consistent with the exact solution when d 2 < 0, does not introduce noticeable errors 
for the values of d 2 that are of concern here. 

The rms deviation over a range of values of d 2 with d\ = is shown in Figure l4~2~l 
(plotted from the exact solution). 

For d 2 = and d\ > the equation can only be solved numerically. The results in 
this case are shown in Figure l4~3l In both cases the deviation from Planckian increases 
approximately linearly with increasing magnitude of the parameters. The simulations 
suggest that for the deviation from Planckian of the CMB to be within the allowed 5 x 
1CT 5 of the peak brightness the diffusion constant d\ must be less than approximately 
7 x 10~ 97 if c?2 = 0, and the drift parameter d 2 must fall within the range — 5 x 10 -96 < 
d 2 < 5 x 10~ 96 if d\ = 0. Converting to SI units this gives the bounds: 

di < 3 x l(T 44 %mV 3 , (4.26) 
-2 x 1(T 43 < d 2 < 2 x 1(T 43 kgm 2 s' 3 . (4.27) 

If both d\ and d 2 are nonzero the picture is somewhat more complicated, but sim- 
ilar bounds result. The results in this general situation are shown in Figure I4.4I The 
allowed region (within the 5e — 005 contour) is seen to be not quite symmetric in the 
drift parameter d%. Placing concrete bounds on d\ and d 2 from Figure |4~41 is difficult, 
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Figure 4.3: The rms deviation of the simulated spectrum from Planckian as a proportion of the 
spectrum peak, varying d\, d 2 = 0. 

since the allowed values of d% depend on the value of c?2 and vice versa - the bounds 
are, however, of the order of magnitude given above. 

4.3.2 The expanding universe 

To determine the bounds on the drift and diffusion parameters, given in Section l4.3.1l 
above, the effect of the expansion of the universe on the CMB was ignored - it was 
assumed that the CMB had remained at a temperature of ~ 2.7 K from the surface 
of last scattering up to today. This is of course not the case. At the surface of last 
scattering the CMB had a temperature of about 3000 K. As the universe expanded 
the individual photons were stretched along with the space, and correspondingly di- 
luted, leaving the 2.7 K spectrum observed today. It turns out that the expansion has 
essentially no effect on the model: the distribution in the expanding universe can be 
deduced easily from the nonexpanding one and the bounds derived from the nonex- 
panding simulation change only slightly. 

Before including the effect of cosmological expansion in the diffusion equation, 
a differential equation that describes the cooling of the CMB needs to be found; the 
following is a standard result in cosmology. Consider the usual continuity equation 

^ = -8^), (4.28) 
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Figure 4.4: The rms deviation of the simulated spectrum from Planckian as a proportion of the 
spectrum peak, varying both d\ and e?2- Values of d\ and within the 5e— 005 contour give a 
spectrum that is Planckian to within 50ppm of the peak brightness. 



for flux pv % . In this case we are interested in the flux in momentum space of the CMB 
photons and thus v has a single component in the energy direction v = dk/dt. This 
gives the continuity equation 



dpt 
dt 







Ok 



dk \ pt dt 



(4.29) 



Under expansion the CMB photons are stretched such that k ~ I /a, where a is the 
scale factor. Thus, dk/dt = —a /a 2 = —ka/a, and the continuity equation becomes 



dpt 
dt 



dk\ a 
adk {ptk) • 



(4.30) 



Consider the distribution 



Po(k,t) 



8tt 



31.2 



a A k 



l o exp 



( ak 
aoTo 



(4.31) 



where ao is the scale factor of the universe today. Photons in the CMB are not only 
stretched, but also diluted by the expansion of the universe. The distribution in Equa- 
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tion 14.311 can be seen as a Planck spectrum where the normalisation depends on the 
scale factor a, to take into account the photon dilution. Equation l4.31l is indeed a solu- 
tion to the expansion equation, Equation l4.30l This confirms, as mentioned earlier, that 
the blackbody nature of the CMB is preserved during the expansion of the universe. 

If the expansion effect is included in the diffusion equation, Equation 14.221 the 
equation becomes 



dpt 
dt 



(di + d 2 )-^- + d x k 




ad. , . 
+ adk M ■ 



(4.32) 



dk 2 



Now consider defining a new variable k = ka/ao and a new density function p(k) 



aoPt(k)/a. Expressed in these variables the spectrum, Equation 14.311 becomes 




(4.33) 



which is constant in time. 



Expressing Equation l4.32l in terms of p and k gives 




RHS 



LHS 



dp 
dt 



(4.34) 



Defining a new time variable t' such that dt' /dt = a/ao reduces this to 



dj_ 

dt' 




(4.35) 



exactly the same form as the non-expanding diffusion equation, Equation l4.22l 

For a matter dominated FRW universe a ~ i 2 / 3 , i.e. a(t) = aot 2 / 3 /t^ 3 , where to is 
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Figure 4.5: The rms deviation of the simulated spectrum from Planckian as a proportion of the 
spectrum peak, varying both d\ and d<i in the expanding universe. Values of d\ and d<i within 
the 5e— 005 contour give a spectrum that is Planckian to within 50ppm of the peak brightness. 

the current value of t. Since dt'/dt = a/ao integration gives 

3t 5 / 3 

t' = -^h + const. (4.36) 

l o 

Since the time since the surface of last scattering (t = 10 60 ) is of the same order as 
the age of the universe to, this effectively gives t' ~ 3/5t. In other words, including 
expansion is essentially the same as running the non-expanding simulations for 3 /5 of 
the time. As such, the bounds found in Section l4.3.1l would not be expected to change 
greatly. For completeness the simulations have been run for the expanding case, and 
the results are shown in Figure [4751 — the order of magnitude of the bounds is indeed 
the same. 
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4.4 The shifting of spectral lines 

The blackbody nature of the cosmic microwave background placed strong bounds 
on the magnitude of the diffusion and drift parameters. It is worthwhile, however, 
looking at other potential astrophysical consequences of the diffusion and drift: are 
observable effects ruled out by our current bounds, or do other sources allow even 
stronger bounds to be derived? Here the effect of drift and diffusion on the spectral 
lines of distant sources will be investigated. 

Diffusion will result in a broadening of spectral lines. The energy dependence of 
the diffusion term in Equation 14 . 1 71 implies higher energy lines will experience greater 
broadening than low energy lines. The difficulty is, however, that spectral lines from 
distant sources will experience line broadening due to a number of other effects. For 
example, the varying velocities of atoms in a gas cloud source will result in thermal 
Doppler broadening of the spectral lines. A calculation of the exact amount of broad- 
ening due to diffusion would not be trivial. Given the very tight bound on the diffu- 
sion parameter d\ < 7 x 10~ 97 , together with broadening due to other effects and the 
inherent spectral resolution (discussed further below), it seems unlikely that broaden- 
ing due to diffusion would be observable. In this section I will instead concentrate on 



The drift term will shift all lines in a spectrum by a fixed amount. Spectral lines 
from astrophysical sources are, of course, already shifted with respect to laboratory 
spectra due to the expansion of the universe. The expansion of the universe results in 
lines shifted by a fixed ratio, the redshift z: 



where X b s is the wavelength we observe, and A e is the wavelength at which the light 
was emitted. Since the drift is independent of wavelength, this effect could be easily 
distinguished from redshifting. In fact, if drift occurs, the calculated redshift for a 
source would depend on which spectral line was used. Assuming that the diffusion 
constant is zero, the drift will shift spectral lines by an energy A A; = d^t, where t is the 
time of travel of the photons. If the line has also been redshifted the observed line will 



drift. 



1 + z 




(4.37) 
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Figure 4.6: The spectrum of quasar SDSSp J083643.85+005453.3 at redshift z = 5.82@ 
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(4.38) 



If drift occurs, but is not taken into account, observations would lead to the calculation 
of an apparent redshift z' that differed from the actual redshift: 



1 + z' 



_ E em n 

Eobs 

, \ E em n 

= (l + Z)— ; r-; — 

V ; E emit + (1 + z)d 2 t 



(4.39) 



As an example, consider the quasar SDSSp J083643.85+005453.3 with redshift z = 
5.82±0.02 (45l, whose spectrum is shown in Figure l4~6l A prominent line in the spectra 
of quasars is the hydrogen Lyman-a line, which has a wavelength of 121.6nm. For a 
quasar at z = 5.82 this line is redshifted to 829. 3nm, or working in energies 



Eemit 
Eredshift 



1.634 x 1(T 18 J, 
2.395 x 10" 19 J. 



(4.40) 
(4.41) 



3 Many thanks to Rick White (Space Telescope Science Institute) and Bob Becker (University of Cali- 
fornia Davis) for providing the data to plot this spectrum. 



; 4.4 The shifting of spectral lines 



91 



To determine how much the Lyman-a line would be shifted by the drift, it is nec- 
essary to know the time it has taken for the light to travel from the quasar. The cosmic 
time-redshift relation 

tu = yft r > ( 4 - 42 ) 

(see for example [46J) gives the approximate age of the universe when photons at a 
redshift of z were emitted. Take the Hubble parameter to be Hq = 70.5kms~ l M-pcr 1 , 
the matter density to be = 0.274 and the current age of the universe to be to = 
13.72 x 10 9 years [47J. The age of the universe when the photons from the quasar were 
emitted is t u = 9.64 x 10 8 years (2.25 x 10 59 in Planck units), and thus they have been 
travelling for t -t u = 1.28 x 10 10 years (2.98 x 10 60 in Planck units). In Section EES the 
drift parameter was constrained to be —5 x 10~ 96 < d,2 < 5 x 1CT 96 . These values give 
a drift of magnitude 

AE = d 2 t 

= 1.49 x 1(T 35 

= 7.30 x 1CT 26 J. (4.43) 



This drift is seven orders of magnitude smaller than the redshifted energy 2.395 x 
10~ 19 J. When measuring the spectra of the quasar, Fan et al. |45] note they have a 
spectral resolution of about 20A. A remeasurement of the spectra |48| gives a much 
better resolution of about 2A. If the measured Lyman-a line is 8293 ± 2A, this corre- 
sponds in energy to approximately 2.395 x 10~ 19 ±lx 1CT 22 J. The drift is thus much 
smaller than even the improved spectral resolution. 

If z = 5.82 is the true redshift of the quasar, the measured redshift from Equa- 
tion 14.391 due to the shift in the Lyman-a line would be imperceptibly different. If 
one were to also consider the Lyman-/? line at 102. 5nm, the difference between the 
redshifts calculated using Lyman-a and Lyman-/? would be of the order 10 -7 , much 
smaller than the stated uncertainty of 0.02. 

The effect of drift would be more apparent for more distant objects. The calcu- 
lations given here for the high redshift quasar at z = 5.82 do however seem to rule 
out any observation of the drift in spectra in the near future. It is certainly clear that 
the bounds placed on the drift parameter through consideration of the CMB are much 
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stronger than any that could currently be obtained from spectral lines. This is of course 
hardly surprising given the combination of the very accurate measurements of the 
CMB and the redshift of z ~ 1000. 
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4.5 Concluding remarks 

The assumption of Lorentz invariance is very powerful: despite the lack of an under- 
lying microscopic model for massless particles, Lorentz invariance allows a concrete 
diffusion model to be derived in the continuum limit. Such models were derived in 
Sections 14.11 and 14.21 in terms of affine time and cosmic time respectively. For mass- 
less particles without internal degrees of freedom, the diffusion model is found to 
have two free parameters - making the phenomenological model, although powerful, 
slightly less powerful than the massive particle case. 

A distribution of massless particles initially peaked at some energy will undergo a 
diffusion in energy at a rate governed by the diffusion parameter d\, and also drift to a 
higher or lower energies depending on the parameter di- The direction of the photon 
is unchanged, and light still travels at the speed of light. 

The values of the two parameters will depend on properties of the, as yet un- 
known, microscopic model. For massive particles the diffusion parameter depended 
on the forgetting time, -ry - for massless particles it is likely there will also be some 
dependence on a non-locality scale, such as that discussed in [39|. The parameters 
could also depend on properties of the wave packet associated with the photon. Fu- 
ture work on scalar field propagation may lead to progress in this area. In the mean- 
time, constraints can be placed on the values of the parameters by comparison with 
observations. 

The CMB, consisting as it does of the 'oldest' photons, provides the perfect labora- 
tory for constraining the drift and diffusion. As shown in Section 14.31 the blackbody 
nature of the CMB constrains d\ and c?2 very strongly. The tightness of the constraints 
is due not only to time the photons have been propagating, but also to the accuracy 
with which the CMB spectrum has been measured. Conceivably, future measure- 
ments could narrow the deviation from blackbody even further and thus result in 
even tighter bounds on d\ and c^. 

The drift and diffusion in energy of massless particles could have other observable 
consequences. One possibility, as discussed in Section l4~4l is the broadening and shift 
of spectral lines. The tight constraints from the CMB do, however, appear to rule out 
the detection of the effect in spectral lines with current technology. 

As mentioned earlier, the work up to this point has considered only particles with 
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no internal degrees of freedom. The ideas here can, however, be applied more gen- 
erally. In particular, as will be seen in the next chapter, the effect of an underlying 
discreteness on photon polarisation can be modelled. 



Chapter 5 



Polarisation 



In the previous chapter the effect of spacetime discreteness on the energy of massless 
particles was discussed, but the possible polarisation of the particles was neglected. 
Here that omission is remedied. 

The study of polarisation is of particular relevance at present - the polarisation of 
the cosmic microwave background (CMB) was first detected as recently as 2002 by 
DASI. It has since been observed by a number of experiments. Current and future 
experiments, such as the recently launched Planck, will hopefully reveal even further 
information. As such, the CMB polarisation is an ideal testing ground for quantum 
gravity phenomenology at the moment. 

If spacetime is discrete it is expected that the polarisation of a photon travelling 
through spacetime will be subject to some fluctuations. Unfortunately there is, again, 
no underlying model for polarisation on a causal set. It is hoped that such a model 
will be developed in the future. In the meantime, the expected continuum behaviour 
of a polarised photon in an underlying discrete spacetime can be found by following 
the stochastic evolution on a manifold of states procedure. 

Before deriving a diffusion equation for polarisation it is necessary to choose, from 
the multitude of possibilities, a suitable way to describe the polarisation. This choice 
is discussed in Section 15.11 In Section 15.21 the affine time polarisation diffusion equa- 
tion is derived. It is found that restricting to linear polarisation is necessary for useful 
concrete phenomenology. The cosmic time linear polarisation diffusion equation is 
derived in Section IB31 and it is found that polarised light will experience both a sup- 
pression in polarisation fraction and a rotation in polarisation angle. 

Section l5~4l discusses possible ways to constrain the free parameters of the model. 
Rough constraints from radio galaxy polarisation data are discussed, as is the effect of 
rotation and suppression on the CMB. 

The majority of the work in this chapter appears in [49 1. 
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5.1 The polarisation state space 

The polarisation of a massless particle with momentum k^ can be described by a com- 
plex four-vector a M , that satisfies k^ = and = 1. This description is not, 
however, gauge invariant: if a M satisfies these conditions then a' 11 = a 1 ' + AA^, for any 
complex number A, will describe the same polarisation state. To avoid this problem, 



the polarisation can be described in terms of a complex two-form P^ u = k^a u 
This P^y satisfies the following Lorentz invariant conditions: 



1. P 



-p 



2. P^P^ = and P* U P^ U = ; 

3. P^k^ = Q; 

4. P^P^ = k v k a . 

Note that it is also the case that any P^ v satisfying these conditions can be expressed 
in terms of some a M in the form above. The antisymmetry condition (1) reduces P^ 
to six complex (12 real) degrees of freedom. The constraints (2) reduce it further to 
10 real degrees of freedom. P^k^ = and P^k^ = give six independent con- 
straints, leaving only four degrees of freedom. Finally, the constraint (4) is only one 
new condition, leaving P^ with three real degrees of freedom. 

Consider the specific case k^ = where = (1, 0, 0, 1). The condition A^a M = 
requires = (a , a 1 , a 2 , —a ), and thus: 
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(5.1) 



where |ai| 2 + |a2| 2 = 1. The overall phase of the complex vector v = (01,02) G C 2 is not 
physically relevant for the description of the polarisation state of a single photon (see 
for example [50 1), leaving a state space with only two real dimensions. This state space 
is the Bloch (or Poincare) sphere, which will be denoted B. The diffusion equations in 
the following sections will be derived in terms of the angles on the Bloch sphere rather 
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than the polarisation two-form P^ u . The resulting equations are of course still Lorentz 
invariant, despite this not necessarily being immediately apparent. For further discus- 
sion of the explicitly Lorentz invariant formulation see Appendix lAl 

The next step is to show that the state space given by v = (a\, CI2) is indeed the 
Bloch sphere. Let m + = 772(1, i) and m_ = -^(1, —%), corresponding to left and right 
circularly polarised states, respectively. The vector v = (01,02) can be expanded in 
terms of {m + , m_}, and an overall phase extracted to leave the coefficient of m + real 
and positive: 
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(5.2) 
(5.3) 
(5.4) 



wherex G [— 7r/4,7r/4] and ip G [0,vr]. Note that this slightly unnatural choice of angles 
is made to match the standard Poincare sphere description of polarisation shown in 
Figure IBTTI (see, for example |51 1). 

The north and south poles (x = ±tt / 4) of the Bloch sphere correspond to 
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right and left circularly polarised states, respectively. States on the equator (x = 0) 
have 
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(5.8) 



and thus the equator consists of the linear polarisation states. The remainder of the 
sphere describes elliptical polarisation. 
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left circular polarization ^ 
Figure 5.1: The Bloch/ Poincare sphere description of polarisation. 

General k^ 

The above calculations showed that for the specific case of k^ = = (1, 0, 0, 1) the po- 
larisation state can be described by the angles \ an d ip on the Bloch sphere. It remains 
to show that for a general k^ the polarisation state P^ can be expressed in terms of 
the Bloch sphere. This can be done by transforming P^ by the Lorentz transforma- 
tion that takes k^ to s M and then simply using the Bloch sphere description for the 
special s M case. In general the Lorentz transformation mapping k^ 1 to s M is not unique. 
A standard Lorentz transformation taking s M to k^ can, however, be defined by (see, for 
example [50]): 

L(k) :=R(k)B z (\k\), (5.9) 

where B z (|k|) is the boost along the z-direction that takes to (|k|, 0, 0, |k|), and R(k) 
is the standard rotation that takes the z-axis into the direction of the unit vector k. If k 
has spherical polar coordinates 9 , 4> then the standard rotation is defined by 



R(k) ■= Rz{(t>)Ry(0) , 



(5.10) 
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i.e. a rotation by angle 8 around the y-axis followed by a rotation by <p around the 
z-axis. The transformation L(/c) _1 takes k^ 1 to s M . Under the same transformation the 
polarisation state, P^ v , gives 

P(k)^ := (LikrXiHkr'VuPp* • (5-11) 

P(k) now has the same form as Equation l5.1[ with components P{k.)^ v = s^a u — a^s u , 
where a\ = P(k)\o and 02 = P(k)2o- Equation 15 . 1 1 1 will henceforth be written P(k) = 
L{k)~ l P for simplicity. 

Thus, a general state (A; M , P^ u ) is specified by coordinates (fc^, VO/ where the 
angles \ and ip are determined by first Lorentz transforming P^ into the form of 
Equation 15. 1[ and then using Equations I5.3l and l5.4l 

5.1.1 Polarisation under Lorentz transformation 

Before moving onto the derivation of the diffusion equation, it is necessary to know 
how the polarisation state (x, ip) changes under Lorentz transformation. The pho- 
ton state (kf 1 , P^u) transforms as usual under a Lorentz transformation A: (k^, P^u) — > 
(&"\ i*V) = ( A ^' 7 ' KK p P v)- For simplicity this will be written (A;', P') = (Ak, AP). 
Under Lorentz transformation (01, 02) — > (a[, a' 2 ) = (P'(k')io, P'{k')2o) and 

P'(k') = L{k')- l P' 

= L(Ak)- 1 AP 

= L{Ak)- l AL{k)P{k) 

= W(A,k)P(k), (5.12) 

where W(A, k) = L(Ak)~ 1 AL{k) is an element of the little group of the Lorentz trans- 
formations that leaves s M fixed. 

The little group is isomorphic to the Euclidean group ISO(2): the set of rotations 
and translations of the plane. Following Weinberg 1501 , an element of the little group 
can be written 

W(9, a, p) = S(a, (3)R(6) , (5.13) 
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where S is the 'translation' 



fi + C ol p -c\ 

a 10 —a 

(3 1-/3 

V C ol P I-CJ 



(5.14) 



a and /? are real numbers and ( = (a 2 + (3 2 )/2. R{6) is the usual rotation 



R${0) 



/l 

cos ( 
sin ( 



(A 
sine 




cos 



(5.15) 



\o o o ly 



To see how the translations 5(a, /3) act on the polarisation consider the infinitesimal 
transformation 

(lap \ 

a 1 —a 
1-/3 
\0 a /3 1 y 

From Equation l5.1l it is clear that AS can be written 



A^(a, P) 



(5.16) 



A5£(a, p) = 8't + s*b v -bl>s v 



(5.17) 



where s M 6^ = and thus W = (6°, a, /3, —6°). Recall P{k)^ v = s^a u — a^s u . Thus 

AS£' A^'P(fc) M v = (< + - b^'s^j (itf + s v 'b v - b v 's v ) 

x (v<v ~ V'V) 



(5.18) 



using the relations s^s^ = 0, s^a^ = and s^b^ = 0. The 'translations' therefore leave 
the polarisation unchanged. 

The rotation R acts as usual, thus under a Lorentz transformation 



(oi, 02) — ► ^2) = (°i cos $ — a 2 sin 5, at sin 5 + a2 cos S) , 



(5.19) 
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where 5 = 6(A, k^). The transformation of the Bloch sphere angles can be determined 
from Equations 15.31 15.41 
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(5.20) 



(5.21) 



thus 



(x, ^) -> (x'> ^0 = (x,^ + 5) , 



(5.22) 



i.e. Lorentz transformations act on the Bloch sphere as rotations around the north- 
south polar axis. 



5.1.2 Stokes parameters 

The Stokes parameters are a common method of describing polarisation, particularly 
in cosmology and astrophysics. They have the advantage of being quantities that can 
be directly measured by experiments. As such, it is useful to understand the relation- 
ship between the Bloch sphere coordinates % , ip and the Stokes parameters I ,Q ,U ,V. 

The first Stokes parameter, /, is the total intensity of the light. The Bloch sphere, 
in fact, only parametrises polarisation of a fixed intensity (or more correctly a fixed 
polarisation intensity Ip where p is the polarisation fraction). The parameters Q and 
U describe linear polarisation. Consider a set of orthogonal axes {x, y} in the plane 
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perpendicular to the direction of propagation of the light, and another set of orthog- 
onal axes {a, b} at 45° in the plane to {x, y}. Q is the difference between the intensity 
transmitted by a polariser that accepts light polarised in the ^-direction, and the in- 
tensity transmitted by one that accepts light polarised in the y-direction. U is given by 
the difference between the intensity transmitted by a polariser that accepts light in the 
a-direction, and the intensity transmitted by one that accepts light in the 6-direction. 
Note that the {x, y} can be chosen freely within the plane perpendicular to the prop- 
agation of the light. In practice there are several standard coordinate conventions, for 
example the International Astronomical Union coordinate convention defines the z- 
axis as pointing along the line of sight towards the observer, the x-axis points north 
and the y-axis, east. The remaining Stokes parameter, V, measures circular polarisa- 
tion. 

The Stokes parameters obey the inequality 

I 2 > Q 2 + U 2 + V 2 . (5.23) 



The polarisation fraction p is defined by 



P =^. (5.24) 



Consider a perfectly polarised beam with fixed intensity normalised to one. The 
relationship between the Bloch sphere angles and the Stokes parameters can be found 
by calculating the density matrix S = vv> from the vector v = (ai, 02). Equating 



1 / I + Q U-iV , 
S = - * , (5.25) 



2 \U + iV I-Q 



(see, e. g. ||52!|), gives 



Q = cos 2ip cos 2x , (5.26) 
U = sin 2ip cos 2 X , (5.27) 
V = sin2x. (5.28) 



More generally, if a collection of photons has a distribution of polarisations given 
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by a probability density p(x, ip) on the Bloch sphere, B, the Stokes parameters are 

Q = P(X, i>) cos 2ip cos 2\d\ dip , (5.29) 
Jb 

U = / p(x,V')sm2V'cos2x(ix#, (5-30) 

v = p(x, fp) sin 2xdxd*p . (5.31) 

Note that a given set of Stokes parameters Q , U , V does not correspond to a unique 
distribution. For example, unpolarised light Q = U = V = could be modelled 
as a uniform distribution of linearly polarised states, or a uniform distribution on 
the two circularly polarised states alone - or any other distribution uniform in ip and 
symmetric in \- 
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5.2 Polarisation diffusion in affine time 

The state space of a photon with polarisation is M. = M 4 x Hq x B. The affine time A 
is again the suitable time parameter for the trajectories. 

Recall that the stochastic process on M can be described by the equation (Equa- 
tion l3"l5) : 

^ = d A {K AB nd B {^)-n A Px ) , (5.32) 

and thus to derive the affine time diffusion equation the Lorentz invariant, symmetric 
positive semi-definite matrix K AB and the Lorentz invariant vector u A need to be 
determined. The microscopic density of states, n, on M. must also be known. 

Unfortunately there are just too many choices when it comes to determining in- 
variant objects on the Bloch sphere. As discussed earlier, Lorentz transformations act 
on the Bloch sphere as polar rotations. Thus, any tensor, vector, or scalar that does not 
depend on is Lorentz invariant. For example, K could be the metric on a sphere, 
or an oblate or prolate spheroid, or indeed any other such strange shape generated 
by rotating a curve around the polar axis. It need not even be invariant under reflec- 
tion in the equatorial plane. The invariant vector, u, can be any linear combination 
°f f(x)d/dx and g{x)d/dij), where / and g are arbitrary functions of \- All these free 
functions mean the model is not much use for phenomenology. Predictive power can 
be regained if the model is restricted to linear polarisation only. For astrophysical and 
cosmological applications this is precisely the area of interest - the CMB, for exam- 
ple, is linearly polarised. The matter is not quite so simple, however, as assuming the 
initial distribution has no circular polarisation. An initially linearly polarised distri- 
bution could become elliptically (or circularly) polarised if K or g{x) are not invariant 
under reflection in the equatorial plane. To allow progress to be made, it will be as- 
sumed from here on that circular polarisation can be neglected. It is hoped that it will 
be possible to determine if this assumption is justified when the microscopic physics 
is better understood. In the meantime, useful results can still be obtained. 

Restricting to linear polarisation, the polarisation state space is just the equator of 
the Bloch sphere: the unit circle, with coordinate tp. In a slight abuse of notation the 
state space with linear polarisation only will also be denoted M. . There is one Lorentz 
invariant vector, d/dijj, on the space of linear polarised states. The Lorentz invariant 
density of states n is constant on the circle. As in Chapter [H it is convenient to work 
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in polar coordinates (k, 9, <p) on the momentum space Hg. The coordinates on the state 
space M will be denoted X A = k, 9, 4>, if)). As shown Section |4~T1 there is the 
invariant vector field k a = (k, 0, 0), and invariant two-tensor k a k b on Hq. The density 
of states on Hq is k sin 9. We also know from Section l4Tl that the spacetime components 
of K AB are zero, and the spacetime component of u A is simply k^. 

Let v\ and V2 be the invariant vectors v A = (0, k, 0, 0, 0) and v A = (0, 0, 0, 0, 1) on 
Ai. The most general symmetric, positive semi-definite invariant two-tensor on the 
full state space M is 



K AB = av A v A + b{v A t 


f + v A 
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(5.33) 



the invariant vector is 



u A = (^,dfc,0,0,e) . 



(5.34) 



Here a > 0, c > 0, b, d, e are real constants. In order for K AB to be positive semi- 
definite the constants must also satisfy ac — b 2 > 0. On A4 the density of states is 
simply n = k sin 9. 

Substituting these forms for K AB and u A into Equation l5.32l above gives the affine 
time linear polarisation diffusion equation 



dpx 



8 A (K AB nd B ( 



-k^+ak 2 
d 2 Px 



n 

d 2 Px 



u A p\ 



+c 



+ 2bk 



0^ + (a-d)k— 
d 2 p\ dp x 



(a + d)p x 



dipdk dip 



(5.35) 



At first glance this equation appears unwieldy, with still too many free parameters 
to be useful for phenomenology. Recall, however, the massless particle affine time 
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diffusion equation, Equation 



dpx ,u d Px , , d f 3 d fpx\\ ,9 



9A dx^ dk \ dk \ k J ) dk 

= ~ k "iB + dik2 ^ + {di ~ d2)k& it ~ {di + d2)px ■ (5 - 36) 



The parameters a and d in Equation 15.351 are just d\ and c?2 renamed. In Section l4~3l 
very tight constraints were placed on the values of d\ and d%. Because d\ and c?2 are 
constrained to be so small, they can be neglected in this discussion of linear polarisa- 
tion. Taking a = d = further requires 6 = 0. Finally the homogeneous polarisation 
affine time equation reduces to 

^ = c ^_ e ^. (5.37) 
d\ dtp 2 dtp 

A linearly polarised photon in an underlying discrete spacetime therefore experiences 
a diffusion, determined by the parameter c, and a drift, given by e, in its angle of 
polarisation. Although Equation 15.371 is Lorentz invariant, the Lorentz invariance is 
not immediately obvious. In Appendix [A] the linear polarisation diffusion equation 
is formulated in an explicitly Lorentz invariant way, using the polarisation two-form 
P^ u . The resulting equation, although equivalent to the equation above, provides less 
immediate physical insight. 
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5.3 Polarisation diffusion in cosmic time 

Again, to compare with observations the affine time polarisation diffusion equation 
must be re-expressed in terms of cosmic time. Following the arguments given in Sec- 
tions [33] and |42] a larger state space incorporating the affine time is defined: M! = 
M. x E. The probability density p\ can be written as a component of the current in M.': 

j\t,x\k a ^,X)= Px . (5.38) 



The t component of the current is 

J t (t,x i ,k a ,ip,X) = kJ x , (5.39) 

where Equation 13.51 has been used. The remaining components of the current, ex- 
pressed in terms of J* are: 

J i (t,x i ,k a ,il},X) = (5.40) 
k 

8 ( ft \ r) ft 

j"(t,x>X,i,,\) = (pr j - ^ + rfJ ' • (5 ' 41) 

= -1% 1*- < 5 - 42 ' 

The unobservable affine time can be integrated over, and a new integrated current 
J defined: 

/•oo 

J t (t,x i ,k a ,^) = / J t (t,x i ,k a ,tp,X)dX, (5.43) 
J o 

J\t,x\k a ^) = ^-J\ (5.44) 
k 

f) / ft \ f) ft 

J\t,x',k\i,) = -at 3 — (^J -6— + <ij', (5.45) 

JUM*) = (J) + «J.. ,5.46) 

The continuity equation on .M' can also be integrated over A to give: 

+ 8^ + dtf + 8 k J k + d^P = . (5.47) 

J x I a=o is zero f° r all i > and J A tends toward zero as A approaches infinity for finite 



J A 
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t. Thus, renaming J* as pt, the polarisation cosmic time equation is 

?£L = - di r - d k J k - d^P 



at 



tfdpt d 2 Pt ,»dpt ^ oh d2 Pt 



k dx l dk 2 dk dkdip 

+ cd^H_2bdp i _edp i 
k dip 2 k dip k dtp 

Neglecting the energy drift and diffusion as discussed in Section \52\ (a = d = =>• b = 
0), and looking at the homogeneous case gives 

d Pt _ Si d 2 Pt 5 2 d Pt 
dt k dip 2 k dip ' 

where, for convenience, the diffusion parameter, c, and drift parameter, e, have been 
renamed 5% and 62 respectively. From this equation we see that the rates of drift and 
diffusion in linear polarisation angle are energy dependent: the higher the energy of 
the photon, the slower the rate of drift and diffusion. This result is particularly inter- 
esting because, although there are other models that predict drift in linear polarisation 
angle, most do not have this particular energy dependence. This will be discussed fur- 
ther in Section l5~4l 

5.3.1 The solution to the diffusion equation 

Equation 15.491 is a standard equation describing diffusion with drift on a circle. Solu- 
tions to this type of equation are well known, see for example II53II . Here I will solve 
the equation using Fourier series. First, note that the diffusion occurs on a half circle 
rather than a full circle: ip G [0, it]. Physically, a polarisation angle ip and a polarisation 
angle ip + tt are equivalent. For simplicity I define a new variable x = ip/ir, < x < 1. 
The distribution p t (t,tp) is a scalar density, rather than a scalar, and thus under this 
rescaling transforms as 

_ , , dip . . 
Pt(t,x) = ■faPt( t ,V) 

= ~Pt(t,iP). (5.50) 

TT 
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Equation l5.49l becomes 
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(5.51) 
(5.52) 



Expand pt in a Fourier series: 

Pt (t, x) 
c n (t) 



E Cn ^ ( 

p t (t,x)e- 2ninx dx. 



n=— 00 
-I 



(5.53) 
(5.54) 



The diffusion equation can thereby be written as an ordinary differential equation 
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u dt 

1 ( _ 62 d &) e - 2 ™ nx dx 
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(5.55) 



Integrating by parts and recalling p t {x = 1) = p t (x = 0) and e n = 1 gives 
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Let the initial distribution be f{x) = pt(0, x), then 
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(5.59) 
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where 



p( x ) = ^2 ex p 



ASin 2 2i5 2 n 
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+ 



A: 



(5.60) 



is the fundamental solution (see for example 153 J). To be consistent with the notation 
that appears in other papers (e.g. |54l |55]|) p can be expressed in a slightly different 
way. Making use of the Jacobi identity 
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p can be written 
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(5.62) 



A; 



47r5it 



S 6XP 



7rn 



(5.63) 



The full solution to Equation l5.49l is thus 
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(5.64) 



(5.65) 
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In the case of zero diffusion, i.e. 5± = 0, the fundamental solution, Equation 15.601 
becomes 



2mn I x 



irk 



(5.66) 



i.e. an initial distribution of 5(x) will become rotated by Szt/nk, corresponding to an 
angle of 82 tj k. For an arbitrary initial distribution f(x), this gives a solution 
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(5.67) 
(5.68) 



Thus, the drift term in the linear polarisation diffusion equation simply rotates the 
plane of polarisation of light by an angle a: ip — > tp' = V 7 + a > where a = S^t/k. In 
terms of the Stokes parameters, V = since there is no circular polarisation and Q and 
U become: 

Q^Q' = Q cos 2a - U sin 2a , (5.69) 
U-fU' = cos 2a + Q sin 2a. (5.70) 

The diffusion term acts as an overall suppression of polarisation. To see this, first 
note that given a collection of photons with a linear polarisation distribution p(tp), the 
Stokes parameters are 

Q = p(ip) cos 2-ipdxl) , (5.71) 
Jo 

f'7T 

U = I p(if>) sin 2ipdip . (5.72) 
Jo 
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Consider a source with an initial distribution pi(tp) = 8(ip). This implies Qi = 
1, Ui = 0, and the light is fully polarised, i.e. the polarisation fraction is p = 1. Suppose 
Pi evolves under some process, such as the diffusion process, to a distribution p(ip), 
the fundamental solution of the process. The Stokes parameters are then 

/*7T 

Q f = / p(ip) cos 2ip dip , (5.73) 

f'7T 

U f = p(tl>)ean2il>dil>. (5.74) 
Jo 

The polarisation fraction is now 

P 2 = Q} + U} 

= p(ip) cos 2^ dip^J + (^J p(ip) sin 2ip dif^J 

< 1 . (5.75) 

Assuming that the process is homogeneous, i.e. 5(ip — ip') — > p(ip — tp'), an arbitrary 
initial distribution p t (0, ip) will evolve to pt{t, ip) = f£ pt(0, ^/)p(ij)—^/)dip' . The Stokes 
parameters are therefore 

/*7T 

Qi = / pt(0,^)cos2^#, (5.76) 
J o 

Ui = f p t (0,^)sin2V>#, (5.77) 

JO 

Q(t) = [ [ p t (0,ip')p(ip-il)')cos2i>dip'di>, (5.78) 
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Jo ■/ o 

Consider the linear polarisation diffusion equation with drift 62 = 0. For zero drift 
the fundamental solution, p(ip), is given by Equation l5.60l with 82 = 0: 
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The Q Stokes parameter is therefore 
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(5.81) 



The second sum vanishes and only the n = ±1 terms in the first sum contribute, 
leaving 

Q(t) = exp (^^j £ Pt (0,ip') cos 2ip'diP' 

= expf— (5.82) 



Similarly, the U Stokes parameter is 
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Here the first sum vanishes and the n = ±1 terms of the second sum give 



U(t) = exp(-^\ £ p t (0,V')sin2W 



expl^^W (5.84) 



The polarisation suppression due to the diffusion is thus 



P 2 



n 2 
Pi 



Q 2 + uf 



(5.85) 



=► p = exp (-^tj . (5.86) 

Both the suppression and rotation are independent of the current polarisation an- 
gle, and thus in the general case of nonzero 5± and 62 a linearly polarised beam in a 
spacetime with an underlying discreteness will experience a rotation in polarisation 
angle and a suppression of its polarisation fraction. 
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5.4 Observational constraints 

The model discussed here is by no means the only mechanism for producing rotation 
of the plane of linear polarisation. Such rotations, when discussed in the context of 
the polarisation of astrophysical or cosmological sources, are usually referred to as 
'cosmological birefringence'. Rotation of polarisation is not merely hypothetical - the 
polarisation from cosmological sources is well known to be subject to rotation, namely 
the Faraday effect. Light travelling through the interstellar medium is subject to ro- 
tation of the plane of polarisation, because the presence of a magnetic field results in 
different refractive indices for left and right circularly polarised modes. The angle of 
rotation, a, is dependent on the wavelength: a a A 2 . 

Proposed methods for generating additional cosmological birefringence include a 
loop quantum gravity motivated modification of the dispersion relation I56U57I . This 
modification leads to circular polarisation modes travelling at different velocities and 
thus can give birefringence. In this case the rotation angle a oc 1/A 2 oc k 2 . Quantum 
gravity motivated birefringence proportional to k 2 is also discussed in [58| and |59|, 
where polarisation data from the Crab Nebula and CMB respectively are used to con- 
strain the model parameter. 

Birefringence also arises in discussions of dark energy, specifically the quintessence 
approach. The coupling of the quintessence scalar field to the pseudoscalar of elec- 
tromagnetism leads to birefringence independent of wavelength l60l l6ll [62)- Vio- 
lating CPT symmetry through the addition of a Chern-Simons term to the Maxwell 
Lagrangian also results in rotation of the plane of polarisation independent of wave- 
length Il63ll64ll65| . Constraints on such polarisation rotation come from radio galaxy 
polarisation data and CMB polarisation. Methods of detecting polarisation rotation 
through radio galaxy data and the CMB will be discussed further below. 

None of the proposed effects discussed above have the a oc 1/k dependence of 
the polarisation diffusion model. One source of a 1/k rotation has been suggested 
previously: according to Prasanna and Mohanty |66| , polarisation rotation could be 
induced by a gravitational wave propagating in the same direction as the photons. 
This effect would only be present where there is a source of strong gravitational ra- 
diation, such as a binary pulsar, and thus it is unlikely to be confused with the effect 
studied here. Note also that most models that predict birefringence violate Lorentz 
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invariance as well as parity. The polarisation diffusion model of this chapter is partic- 
ularly interesting because it preserves Lorentz invariance. 

Before discussing the bounds on the model parameters, it should be noted that 
special care needs to be taken when searching for changes in polarisation. Unless the 
method of generation of polarisation in a source is completely understood, natural 
variations in polarisation at the source could be mistaken for effects occurring over 
the travel time of the light to us. 

5.4.1 Constraining the drift parameter with radio galaxy data 

Radio galaxies emit linearly polarised light and often have an elongated structure 
defining a natural axis against which their polarisation angle can be measured. Af- 
ter subtracting the effect of Faraday rotation it is found that the polarisation angle for 
radio galaxies is highly correlated with the observed galaxy axis (radio axis), peak- 
ing at 90° to the radio axis (see, for example, Il67ll64l I68H and Figure I5T21) . The linear 
polarisation of radio galaxies is due to the synchrotron process. Although the cause 
of the correlation between polarisation angle and radio axis does not seem to be fully 
understood it is likely there is a physical basis for the alignment of the magnetic field. 
The correlation holds over a range of redshifts and thus provides a method of con- 
straining a birefringence effect. If it is assumed that the light has a polarisation of 90° 
at the source, the difference from 90° in the measured polarisation angle constrains 
the amount of rotation, and in the case of the model discussed here, the magnitude of 
the drift parameter. Data also exists for polarisation fraction for radio galaxies. Un- 
fortunately there seems to be less understanding of what the polarisation fraction is at 
the source and thus the change in polarisation fraction (and the diffusion parameter) 
cannot be constrained. 

In this section I will use data from one particular galaxy to place an order of mag- 
nitude bound on the drift parameter. In II68II the galaxy MRC 2025-218 was stud- 
ied. Cimatti et al. deduced a polarisation angle (with respect to the galaxy axis) of 
9 = 93 ± 8° and a polarisation fraction of p = 8.3 ± 2.3%. They concluded that there is 
"strong evidence that the plane of polarisation is not rotated by more than 10 degrees 
when the radiation travels from z = 2.63 to us". Recall from Section 1531 that the angle 
of rotation is related to the drift parameter by a = 52t/k. The time of travel t can be 
calculated from the redshift, z = 2.63, of the galaxy. As mentioned in Section l4~4l the 
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Figure 5.2: The correlation between radio galaxy axis and polarisation angle. Data from 71 
galaxies with redshift z > 0.3 is shown here. This data is taken from l64ll67l . 



cosmic time - redshift relation, 



(5.87) 



gives the approximate age of the universe when photons at a redshift of z were emit- 
ted. Taking the Hubble parameter to be Hq = 70.5fcms _1 Mpc~ 1 , the matter density 
to be Qm = 0.274, and the current age of the universe to be to = 13.72 x 10 9 years 
(see [47]), gives t u = 2.55 x 10 9 years, and thus the photons have been travelling for a 
cosmic time oft = 11.17 x 10 9 years. 

Cimatti et al. observed the galaxy MRC 1015-218 in the R band, corresponding to a 
rest frame spectral region (i.e. non-redshifted wavelength) of AA res t = 1600 to 2100A. 
In the calculation below, I will assume the photons were emitted with a wavelength 
of A = 1850A. As the polarisation rotation angle, a, depends on the wavelength, some 
error is introduced by taking A to be a fixed value within the bandwidth. The variation 
of a over this wavelength range is not expected to be significant however, and should 
not affect the order of magnitude bound obtained. 

The polarisation diffusion model assumes Minkowksi spacetime, and thus to cor- 
rectly bound the drift parameter the change in energy of the photons from the galaxy 
due to the expansion of the universe must be taken into account. The inclusion of 
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expansion into the model will be discussed more generally in Section l5.4.2l here I will 
only note that the wavelength of photons in an expanding universe is given by 

A(t) = A e ^ , (5.88) 

where A e is the wavelength the photons were emitted at, a(t) is the scale factor, and 
a e is the scale factor at the time when the photons were emitted. Assuming a matter 
dominated FRW universe, the scale factor goes as a ~ t 2 / 3 , i.e. a(t) = t 2 / 3 /t 2 /3 where 
to is the current age of the universe and the scale factor today is taken to be ao = 1. 
The total angle of rotation of the polarization is thus 



Aa 




(5.89) 



The drift parameter can be constrained to be 



So < 7 x 1(T 90 



2 

37r 2 -3 



< 2 x 10~ at kgm z a~ a . (5.90) 

5.4.2 Polarisation diffusion and drift in the cosmic microwave background 

The effect of polarisation rotation on CMB data has been well studied. Suppression 
of polarisation, however, seems to have been little considered. Before discussing the 
CMB we need to fully understand how the expansion of the universe changes the 
diffusion equation. In fact, an equation that includes expansion can be determined in 
a similar manner to that discussed in Section 1431 
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Including the expansion of the universe 

The redshifting of the photons is given by 

k(t) = k e ^f- (5.91) 
a(t) 

where k e is the emitted frequency, a(t) is the scale factor, and a e is the scale factor at the 
time the photons were emitted. Note that the expansion of the universe also dilutes 
the photons, but this dilution does not change the diffusion equation. As mentioned in 
the previous section, if a matter dominated FRW universe is assumed, a(t) = t 2 ^ 3 /t^ 3 , 
where to is the current age of the universe and the scale factor today is taken to be 
do = 1. If the form for k is substituted into the diffusion equation, Equation 15.491 it 
becomes 



dp t 5ia(t) d 2 p t 5 2 a(t) dp 



dt k e a e dip 2 k e a e dip 

3 ( Sx d 2 Pt 6 2 d Pt 



t 2 / 3 \k e a e da 2 k e a e dtp 



(5.92) 



A new time variable t' can be defined by 



dt' 

— = a(t), (5.93) 



and thus 



dp t t 2 / 3 dp t 



dt 1 t 2 / 3 dt ' 

dp t dp t dt! 
dt ~ dt' dt 
t^dpt 

,2/3 dt' ' 

dp t 5i d 2 p t S 2 dp t 



(5.94) 



dt' k e a e dip 2 k e a e dip 



(5.95) 
(5.96) 



If we further note k e /ko = ao/a e , where ko is the energy observed today, the equation 
becomes simply 



dp t Si d 2 p t 5 2 dp 



dt' ko dip 2 ko dip 



(5.97) 
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In other words, the change in energy due to expansion can be compensated for by fix- 
ing the energy at the observed energy and simply transforming to a new time variable 

3t 5 / 3 

t' = --^r + const. (5.98) 

For analysis of the CMB this is approximately equivalent simply to evolving the equa- 
tion for 3/5 the time since t is so close to to- 

The polarisation of the CMB 

The cosmic microwave background is expected to be partially polarised due to Thom- 
son scattering off free electrons at the surface of last scattering. Thomson scattering 
can produce only linear polarisation (and only do so when there is a quadrupole 
anisotropy in the incident radiation) and thus the CMB is expected to have no circular 
polarisation. The CMB polarisation has now been detected by multiple experiments, 
including WMAP EH, BICEP (701, and QUaD £fl. 

The polarisation of the CMB is commonly discussed not in terms of Stokes param- 
eters, or angles on the Bloch/Poincare sphere, but rather in terms of E and B modes 
(see, e. g. II72II ). Instead of describing the polarisation of a single photon, E and B 
modes describe how the polarisation varies over a small area. The names 'E' and 'B' 
mode are drawn from an analogy with electric and magnetic fields: the polarisation 
pattern on the sky (direction and intensity) is decomposed into a curl-free component 
(E mode) and a divergence free component (B mode). The exact method by which 
these components are calculated is not important for the current discussion. As men- 
tioned above, a quadrupole anisotropy in the radiation at last scattering is necessary 
for there to be polarisation in the CMB. This anisotropy can arise from scalar, vec- 
tor, or tensor perturbations. Scalar perturbations, due to density fluctuations in the 
plasma, give rise to E mode polarisation; vector perturbations, due to vorticity in the 
plasma, are expected to be negligible; tensor perturbations can be caused by gravita- 
tional waves, and would give rise to B mode polarisation. 

Analysis of the CMB is usually in terms of correlations between E mode, B mode, 
and temperature T. In all, there are six cross-correlation spectra that can be examined: 
TT, TE, TB, EE, EB, BB. In the absence of a parity violating effect, the correlators 
TB and EB are expected to vanish. The rotation of the polarisation angle can thus 
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be detected through the presence of a nonzero TB and EB. More specifically, for an 
angle of rotation Aa the exact change in the correlators can be written as (see, for 
example E2UZ2UE3): 

C'f B = Cf s sin 2Aa, (5.99) 

cf B = \{cr-cr)^^, (5.100) 

C'f E = Cf E cos2Aa, (5.101) 
C'f E = Cf £ cos 2 2Aa + Cf B sin 2 2Aa , (5.102) 
C'f B = Q BB cos 2 2AQ + Cf £; sin 2 2Aa, (5.103) 

where C\ are the spectra when there is no rotation. Note that the TT correlation is 
unchanged by rotation. For the polarisation diffusion and drift model discussed here 
the angle Aa is dependent on the energy, and thus the expected spectra above will 
depend on the frequency band one is working in. 

Example correlation spectra can be produced using the Code for Anisotropies in 
the Microwave Background (CAMB) ] Figure E3 shows TT, TE, EE and BB spectra 



1 Available at http:/ / camb.info/. 
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(b) EB. 

Figure 5.4: CMB TB and EB correlation spectra when polarisation is rotated by an arbitrary 
angle Aa = tt/8. 

for standard parameters, without any polarisation rotation (recall TB and EB van- 
ish). Figure [5~4l shows the TB and EB spectra that would result if a rotation of po- 
larisation occurs; an arbitrary angle Aa = tt/8 has been chosen to illustrate the effect. 
The changes in TE, EE, and BB spectra under this rotation are shown in Figure 1531 
Note that there is some difficulty determining what values of Aa are allowed by cur- 
rent data. For example, Xia et al. H74l give differing (and seemingly contradictory at 
68% C.L.) values of Aa = -2.62 ± 0.87°, Aa = 0.59 ± 0.42°, and Aa = 0.09 ± 0.36° 
depending on which subset of CMB experiments they analyse. Note that as they as- 
sume Aa is frequency independent these values cannot be used to constrain the drift 
parameter 62- Recall that a = Syt' /k when expansion of the universe is included. Tak- 
ing t' = 3/5 x 10 60 , and choosing one of the WMAP observation frequencies, 40GHz, 
shows that the (probably too large) angle Aa = tt/8 in fact corresponds to 62 ~ 10 -93 , 
already several orders of magnitude smaller than the bound placed on 82 with radio 
galaxy data in Section 15.4.1 1 
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Figure 5.5: CMB TE, EE, and BB correlation spectra when polarisation is rotated by an ar- 
bitrary angle Aa = ir/8 (dashed lines) compared to the spectra in the unrotated case (solid 
lines). 
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The discussion above has all been concerned with the effect of the rotation of the 
plane of polarisation on CMB data. Effects such as the suppression of polarisation 
predicted by the polarisation diffusion equation do not seem to have been discussed 
much in the literature. In terms of the correlation spectra above, the diffusion will sim- 
ply reduce the amplitude of each spectra by some fixed, energy dependent, amount. 
To properly constrain the drift and diffusion parameters, CMB data from a single fre- 
quency band must be compared with the corresponding predicted spectra generated 
by CAMB and modified for drift and diffusion. Comparing over a range of drift and 
diffusion parameters will allow the values that best fit current CMB observations to be 
found. Unfortunately this process requires modifying the Markov Chain Monte Carlo 
code, CosmoMC, and is outside the scope of this thesis. This work has recently been 
undertaken by Carlo Contaldi and the results appear in [49|. 
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5.5 Concluding remarks 

When studying massless particles with polarisation in a discrete spacetime, the as- 
sumption of Lorentz invariance turns out to be insufficient to completely constrain 
the phenomenological model. To make concrete predictions it is necessary to restrict 
attention to linear polarisation - fortunately it is precisely linear polarisation that is of 
interest in astrophysics and cosmology. 

As yet there is no underlying microscopic model for polarised photons on a causal 
set (or even, as discussed in Chapter [4] unpolarised photons). When particles in causal 
set theory are better understood it is hoped that useful models including circular po- 
larisation will be able to be developed and the validity of the results obtained in the 
current restriction to linear polarisation checked. 

For linear polarisation a polarisation diffusion equation was derived in Sections l5.2l 
and 15.31 in terms of affine and cosmic time respectively. The tight constraints on en- 
ergy diffusion and drift obtained in Chapter [4] allow the polarisation diffusion equa- 
tion to be simplified to a diffusion equation in terms of polarisation angle with two 
free parameters. A distribution of linearly polarised photons is found to experience a 
suppression in polarisation determined by the diffusion parameter Si and a rotation 
of polarisation determined by the drift parameter 82 ■ 

The suppression of polarisation over the travel time of photons from distant sources 
to us does not seem to have been much considered. Unfortunately it is, for the most 
part, a difficult effect to constrain since the polarisation fraction at the source must 
first be well understood. 

There are a number of different proposed models that would generate polarisation 
rotation over the travel time of photons. The model proposed here is distinct, how- 
ever, in the l/k dependence of its rotation and in that it preserves Lorentz invariance. 
The polarisation rotation is easily constrained by existing data. A natural correlation 
between the polarisation angle and the axis of radio galaxies allowed a constraint of 
62 < 2 x 1CT 37 kgm 2 s~ 3 to be placed in Section l5T4l The photons from the CMB, having 
had more time to experience the underlying discreteness, impose a tighter constraint 
on the rotation and also constrain the suppression. These constraints are discussed 
in 11491. 



Chapter 6 



Conclusions and future 
directions 



Lorentz invariance is not optional in the causal set approach to quantum gravity. It is 
an essential feature of the theory and proves a very powerful constraint in the devel- 
opment of causal set phenomenology. A surprising amount can be deduced beginning 
with only the simple assumptions that spacetime is discrete and Lorentz invariant. 

The behaviour of massive particles in a discrete spacetime was shown to be de- 
scribed by a diffusion equation with one free parameter in the continuum limit. The 
derivation of this equation did not rely on any underlying particle model, but sim- 
ply on the assumptions of discreteness and Poincare invariance. It could be argued 
that the derivation does not guarantee a finite diffusion constant and thus a truly ob- 
servable effect. To rebut this objection it was shown analytically that the continuum 
diffusion parameter was necessarily finite for the swerves model of particle propaga- 
tion on a causal set. Numerical simulations of the swerves model and two 'intrinsic' 
models gave results that clearly demonstrated diffusion behaviour, and in the case of 
the swerves model it was possible to determine the relationship between the diffusion 
parameter and the underlying model parameters. 

In the causal set approach, spacetime discreteness is expected to be of the order 
of the Planck scale. The constraints of computer memory limit causal set sizes simu- 
lated on a standard computer to around 2 17 elements. We could question the point of 
even attempting to investigate large scale behaviour in causal sets that correspond to 
such absurdly small regions of spacetime (even a volume of 1mm 3 x Is would con- 
tain on the order of 10 137 elements). Yet remarkably, the simulations do agree with 
the expected continuum behaviour. The results suggest that causal set simulations 
will have an important role to play in the further development of the theory and its 
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phenomenology. 

Lorentz invariance also allowed a useful phenomenological model for massless 
particles in discrete spacetime to be developed. Massless particles were shown to ex- 
perience a drift and diffusion in energy in the continuum limit. Crucially, all massless 
particles still travel at the speed of light in this model. The accurate data available on 
the cosmic microwave background allowed both the drift and diffusion parameters to 
be tightly constrained. Whether the constraints are so tight as to exclude all possibility 
of future observations of the effect remains to be seen. 

The common word 'particle' erroneously suggests massless and massive particles 
are variations of the same fundamental object. The diffusion experienced by massless 
particles is not simply the massless limit of the massive particle diffusion, suggesting 
that massive and massless particles will appear in quite different forms in any future 
complete causal set theory. How exactly will particles fit into the causal set approach? 
The massive particle models considered here were classical, point particle models and, 
as yet, we have no candidate models at all for massless particles. To produce any 
underlying model for massless particles, it seems likely quantum behaviour will have 
to be considered. 

Placing particles 'on' a fixed causal set seems highly unsatisfactory. The causal 
set should be a dynamic structure, growing in number of elements. On large scales 
the theory must reproduce the relationship between mass and curvature, suggesting 
that the presence of particles will influence the growth of the causal set. The need 
to reproduce quantum effects suggests assigning amplitudes to each possible step in 
the causal set growth. It would be nice if we could just devise a theory now that 
satisfies all our requirements, but the problem is simply too big. Progress is made by 
looking at what can be learned from phenomenology, by developing quantum field 
theory on causal sets, by looking at the structure of curves spacetime in causal sets, 
and so on. Steps in the direction of a complete theory will hopefully lead to further 
phenomenology and, ultimately, falsifiable predictions. 

The power of the Lorentz invariance assumption has its limitations. When con- 
sidering the polarisation of massless particles, Lorentz invariance proved insufficient 
to constrain the model. Identifying the assumptions necessary to constrain the model 
further will have to wait for the next step in understanding particles in causal set 
theory. In the meantime, restricting to linear polarisation has allowed progress to be 
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made. In a discrete spacetime linearly polarised photons were found to experience 
rotation in polarisation angle and suppression in polarisation fraction. Rotation of po- 
larisation angle is an effect predicted by a number of other approaches to quantum 
gravity but the particular energy dependence of this model is unique. Again, cosmo- 
logical data allowed the free parameters of the model to be constrained. 

It is important to remember that the results discussed here do not, for the most 
part, call on specific features of causal set theory. They rely on the assumptions of 
Lorentz invariance and discreteness and may, hopefully, be generalisable to other 
quantum gravity theories. A lot has been achieved looking at the large scale behaviour 
of particles in discrete spacetime, despite the incomplete nature of the theory. Where 
next can we look for signatures of spacetime discreteness? 

The importance of phenomenology in the search for a theory of quantum gravity 
should not be underestimated. It is crucial in the development of the theories and, 
more importantly, reminds us to look out at the universe once in a while and remem- 
ber why it is we seek such a theory to begin with. 



Appendix A 



Explicitly Lorentz invariant 
polarisation diffusion 



The state space for a photon with polarisation is Ai = M 4 x Iq x P, where P is the 
polarisation state space. For an explicitly Lorentz invariant formulation of the polari- 
sation diffusion equation, P is the polarisation state space with coordinates P^ u , where 
Pfiv = kfj,a u — CLpky is the complex two-form satisfying the conditions set out in Sec- 
tion 15.11 If we wish to discuss only linear polarisation, P^ v can be taken to be a real 
two-form satisfying the conditions 

IP — —P ■ 

2. P^P^ = 0; 

3. P^ = 0; 

4. P^P^ = k v k a . 

Again, any P^ v satisfying these conditions can be expressed as P^ u = k^a v — a^ky, 
where k^d 1 = and a^a^ = 1. a M is now a real four-vector. Condition (1) above 
provides six constraint equations, condition (2) a further one constraint. Although 
condition (3) is four equations, it gives only three new independent constraints. Fi- 
nally condition (4), although 16 equations, gives only one new constraint. The state 
space P therefore has only one degree of freedom. 

A.l Invariant vectors 

Before writing down the diffusion equation, K AB , u A , and the density of states n must 
be known. Consider the submanifold Hi] x P. There are two Lorentz invariant vector 
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fields on the submanifold: vf = (k fl ,P flu ) and v£ = (0,*P^), where *P flu is the 
Hodge dual defined in the usual way as = ^-eP a ^ u P pa . To see this, first consider 
the more general vector fields = (aik^, /3iP fiV ) and = (a^fe^, fh * P^ u )- These 
vectors are clearly Lorentz invariant, as is the submanifold Hq x P. It remains to show 
that the vectors are tangent to the submanifold. 

The normals to the submanifold are derived from the constraint equations in the 
usual way: 



=0 n\ 



d{k v k v ) d(k v k v ) 



P^Pnv = 



[XV — kj —r 71 A 



A \ dk>* ' dPP° 

(2^,0), (A.l) 

2 _ fd(p^p pa ) djp^p^) 



Q k H ' Qppa 

(0,2P pa ), (A.2) 
d(P ua k u ) d{P va k u ) 



" k " ° s dW ' 8PP° 

= (P„ a ,6 a ak p ) , (A3) 

P^Pfw = k v k a => n^ A = (j^; {P^Pfw — k u k a ) , Qp a/3 {P^P^a ~ k v k a )^ 

= {-g. a k v - %k a , b^P aa + g aP P a v ) . (A.4) 

The vector w\ is clearly orthogonal to the first three normals: 

n A w i = 2aik p kP 

= 0, (A.5) 
n\w£ = 1$ x P pry PP° 

= 0, (A.6) 
n'Xwf = a 1 P lx a k^ + p 1 6 a a k p PP' T 

= 0. (A.7) 

For the last normal 



n 



» A wf = (-g aa k u - 8 v a K) a,k a + (8 v p P aa + g a pP a v ) fti**" 
= -ai (k"k a + k v k a ) + 0x {P aa P au + P a v P a c ) 
= 2 (- ai k v k a + hP aa P av ) . (A.8) 
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This equals zero only if ol\ = f}\. Thus the vector v A = (/c^,P pCT ) is normal to the 
surface. 

w 2 is also obviously orthogonal to the first three normals if it is first noted that 
p^ * pp* = o and *P^ u k fl = if both P^P^ = and P^k^ = 0. These results can 
be shown simply if P pu is expressed as P^ u = k^a v - a p k u . 

n A w 2 = 2a 2 k p k' x 

= 0, (A.9) 

n> 2 A = 2p 2 P p<T *PP° 

= 0, (A.10) 

n a A w A = a 2 P^k» + f3 2 5 a a k p *P pa 

= 0. (A.11) 



The remaining normal gives 



<A W 2 = {~9acM v - 5 v a K) a 2 k a + (5%P aa + g^P a u ) fo * P 



= -a 2 {k v k a + k v k a ) + f3 2 (P aa * P av + P a v * P a a ) 

= 2 {-a 2 k v k a + foP au * P av ) . (A.12) 

It turns out that P aa * P au = 0, this is easily checked by choosing a specific frame, 
say k^ = (k, 0,0, k). Therefore w 2 is only tangent to the surface if a 2 = 0. Thus 
v 2 = (0, *P^ V ) is an invariant vector. 

A.2 The diffusion equation 

Given the above invariant vectors v\ and v 2 , the most general Lorentz invariant, sym- 
metric, positive semi-definite matrix K AB on the submanifold Hq x P is 

K AB = av A v B + b (v A v B + v£v?) + cv A v B , (A.13) 

the general invariant vector is 



u A = dv A + ev 2 , 



(A.14) 
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where a, b, c, d, e are constants and ac — b 2 > 0. The spacetime components of K AB 
and u A on the full state space Ai are unchanged from the zero spin massless particle 
case: K pA = 0, u p = jfe**. 

The diffusion equation can thus be written down from Equation l3.15l 

— = d A (K nd B [-)-u Px 



d 

+ 



+ 



d 



akW^- ( —) + (ak p P pa + bkP * P pa ) n-^— (—)- dk" Px 
dk v \n J y ' dPP° \ n 1 

(ak p P pr7 + bk p * P pa ) n-£- ( ?± 



+ {aP pC7 P^ + b (P pa * P^ s + *P^PT" 5 ) + c * P pa * P^ n-^p^ (— 



{dP pa + e * P pa ) p x 



(A.15) 



Although it explicitly demonstrates the Lorentz invariance of the process, this 
equation is hardly user-friendly. To understand the physical consequences of the dif- 
fusion, the equation derived in Section 15721 using the Bloch sphere as the polarisation 
state space is far more useful. It is important to confirm however, that these equa- 
tions are in fact equivalent. First consider polar coordinates onHg: {k,0,<p}. In these 
coordinates the vector k p has components (k, 0, 0). As mentioned earlier, the polarisa- 
tion state space P has only one degree of freedom - specifically, it can be shown that 
the constraints on P pv reduce the polarisation state space to a circle. The density of 
states, n, on this circle is simply a constant and thus the density of states for the full 
state space M. is n = k sin 6. The submanifold Hjj x P is thus the momentum space, a 
null cone, with a copy of the polarisation space, a circle, at each point. The invariant 
vector v A = (k p , P^ u ) is simply the vector at a point on the cone that points along 
the cone, away from the tip. The invariant vector v A = (0, *P pu ) is the tangent vector 
to the polarisation circle. Denoting the coordinate on the circle ip, v A and v A become 
v A = (k, 0, 0, 0) and v A = (0, 0, 0, 1) in coordinates {k, 9, (f>, ip}. Thus, this formalism 
is exactly the same as that used in Section 15721 
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